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We consider Af = 2 supersymmetric gauge theories perturbed by tree level superpo- 
tential terms near isolated singular points in the Coulomb moduli space. We identify the 
Seiberg-Witten curve at these points with polynomial equations used to construct what 
Grothendieck called "dessins d'enfants" or "children's drawings" on the Riemann sphere. 
From a mathematical point of view, the dessins are important because the absolute Ga- 
lois group Gal(Q/Q) acts faithfully on them. We argue that the relation between the 
dessins and Seiberg-Witten theory is useful because gauge theory criteria used to distin- 
guish branches of M = 1 vacua can lead to mathematical invariants that help to distinguish 
dessins belonging to different Galois orbits. For instance, we show that the confinement 
index defined in hep-th/0301006| is a Galois invariant. We further make some conjectures 



on the relation between Grothendieck's programme of classifying dessins into Galois orbits 
and the physics problem of classifying phases of M = 1 gauge theories. 
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1. Introduction 

The interest in Af = 2 supersymmetric gauge theories is especially due to the seminal 
work of Seiberg and Witten in the mid nineties |l|,[2[]. They showed that the low energy 
action and infrared dynamics of the gauge theory on the Coulomb branch can be completely 
solved and that all the relevant information about the low energy theory is encoded in a 
hyperelliptic curve and in an associated meromorphic differential. This work led to a 
tremendous amount of progress in the understanding of the physical aspects of Af = 2 
gauge theories, including the vacuum structure of related Af = 1 theories. At the same 
time, there have also been fascinating connections between Seiberg- Witten theory and 
mathematics, especially to the Donaldson theory of four manifolds @. 

In this article, we unearth a new connection between a particular class of Seiberg- 
Witten curves and what Grothendieck called "dessins d'enfants" or "children's drawings" . 
We will refer to these simply as "dessins" in what follows. We will formally define a dessin 
later on, but for the moment by a dessin we simply mean a connected graph on a two 
dimensional surface, with vertices of two kinds - say black and white - that alternate along 
the graph. 

The original reason for studying such drawings in mathematics was that there is a 
natural action of the absolute Galois group Gal(Q/Q) on them. Moreover, the action is 
faithful, i.e. there is no group element, other than the identity, that leaves invariant all 
dessins. The absolute Galois group is one of the central objects of interest in mathematics, 
especially in number theory. Surprisingly, this object has already made its appearance 
in the physics literature in the context of rational conformal field theory due to work by 
Moore and Seiberg j|,|5],[| . It has been known that the solutions of the Moore-Seiberg 
equations lead to a projective representation of the so called Teichmuller tower. As noted 
by Grothendieck in 0, the absolute Galois group also acts on this tower. Amusingly 
enough, both the Teichmuller tower and the dessins were introduced by Grothendieck in 
the same letter J7|. 

In the rest of the introduction we will summarize recent progress made in ||[J| in 
understanding the vacuum structure of Af = 1 gauge theories obtained by deforming 
an Af = 2 theory by a tree level superpotential. We will see that these gauge theory 
techniques and the existing results could have implications for the study of the action of 
the Galois group on the dessins and, conversely, we would like to argue that there is a 
wealth of information on the mathematical side that could potentially lead to an improved 
understanding of the phases of Af = 1 gauge theories. 
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More generally, we conjecture that there is an intimate relation between Grothendieck's 
programme of classifying dessins into Galois orbits and the physics problem of classifying 
certain special phases of M = 1 vacua. The precise form of the conjecture is given in 
Section 4.4. The meaning of the different physical and mathematical elements that go 
into this conjecture will be explained in the remaining part of Section 1 and in Section 2, 
respectively. 

This paper is organized as follows: In Section 2, we discuss the dessins and related 
mathematical material. This includes Belyi maps, the action of Gal(Q/Q) on dessins 
via Belyi maps and the correspondence with Seiberg-Witten curves. In Section 3, we 
review the concept of invariants (order parameters) which are used to distinguish orbits 
of dessins under the action of Gal(Q/Q) (Af = 1 gauge theory phases). In Section 4, 
we give examples of cross fertilization that arise from our identification of the physics and 
mathematics programmes. In particular, we explain how the confinement index introduced 
in H leads to a Galois invariant. In Sections 5 and 6, we give some illustrative examples in 
full detail. In Section 7, we conclude with many open questions and directions for future 
work. In Appendix A, we give a very basic introduction to algebraic concepts in Galois 
theory with the goal of understanding the definition of Gal(Q/Q) and its action. In the 
same appendix, a small glossary of terms used in the text is given. Finally, in Appendices 
B and C we expand on some interesting issues discussed in the rest of the paper. 

1.1. Physics Preliminaries 

In the Seiberg-Witten solution of Af = 2 gauge theories it is natural to study the loci 
in the moduli space where the SW curve develops singularities. For a pure U(N) gauge 
theory - the most well studied case - the curve has the form 



where Pn( z ) = (det(zl — $)) and $ is the adjoint scalar in the vector multiplet. (We 
adopt the usual notation, with the subscripts denoting the degree of the polynomials.) A 
well studied || factorization is 



which corresponds to a family of degenerate curves. Let us briefly discuss the physical 
information that is encoded in this kind of polynomial equation. 
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The Af = 2 Coulomb moduli space of the gauge theory is parametrized by iV param- 
eters Uk = -^(Tr $ fc ) , constructed from the adjoint scalar $. The problem studied in the 
physics literature is to find and classify the Af = 1 supersymmetric vacua obtained by 
perturbing the Af = 2 theory by a tree level superpotential 

n+l 
fc=l 

Once the tree level potential is introduced, all points in the Coulomb moduli space are lifted 
except those for which N — n mutually local magnetic monopoles become massless. The 
superpotential triggers the condensation of monopoles and the magnetic Higgs mechanism 
leads to confinement of the electric charges. The points at which this occurs are precisely 
those that solve the factorization problem (|1 . 2|) . At low energies at these points, out of 
the oris; inal U(1) N only a U(l) n sub group remains unbroken and its coupling constants 
are given by the reduced Seiberg-Witten curve 

y 2 = F 2n (z). (1.4) 



From a simple counting of parameters, we see that ( |1.2| ) defines an n-dimensional 
subspace of the moduli space. Plugging the u^s as functions of these n parameters into 
the superpotential Wtree , one gets an effective superpotential 



n+l 



W eS = S ^2g k u k (1.5) 



fc=i 

on the moduli space. Thus, given the parameters g k , one can vary with respect to the 
coordinates on the moduli space and obtain all the u^s as functions of the g^s and A. 

In other words, the form of the superpotential picks out specific points in the Af = 2 
moduli space which correspond to Af = 1 vacua. What was shown in JTD| , [TT[] is that this 



extremization problem can be rephrased as the problem of factorizing the Seiberg-Witten 
curve in the following manner: 

P 2 N {z) - 4A 2Ar = -^{{WLe(z)f + U-i(z)) H 2 N _ n (z) . (1.6) 

9 n+l 

Let us denote the critical points of W tre e by a^, i.e., W[ ree (z) = Y\7=i 9n+i(z ~ a i)- We set 
<7n+l = 1 in what follows. 
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It is interesting to ask what the "moduli space" of M = 1 vacua is, as the parameters 
<7fc's are varied. The motivation for such a question was explained in ||: semiclassically, 
as A — > 0, the gauge group U(N) can be broken to niLi U(Ni) with = N by 

choosing $ to be a diagonal matrix with N{ entries equal to It is therefore natural to 
ask whether it is possible quantum mechanically to interpolate between vacua that have 
the same n but different iVj's. 

This was answered in Q where it was shown that, for example, vacua with one U(l) 
factor can be smoothly connected to vacua with any allowed values of iVj's. All these 
classical limits are different corners of a single connected subspace of the M = 2 moduli 
space, which was referred to as an M = 1 branch. However, in || it was also discovered that 
there were other branches distinguished by order parameters such as the expectation values 
of Wilson loops. Branches meet at points where extra massless mutually local monopoles 
appear. At these points other branches also emanate which have the same dimension but 
where some of the N^s are zero. 




Figure 1: Two J\f = 1 branches of the U{4) gauge theory obtained by deforming the Af = 2 
theory by a cubic superpotential. Only those branches with two U(l) 's in the infrared are shown. 
The two branches meet at a point where one more monopole becomes massless. 

The structure of these branches and how they meet can be quite intricate ||. We 
use the example of a cubic superpotential in a U(4) gauge theory to represent a region 
close to one of the points with three massless monopoles in Figure 1. At a generic point 
the low energy group is U(l) 2 . At the special point where the two branches meet, a third 
branch (not drawn) emanates which consists of vacua with a single U(l) as the low energy 
group. Along each of the depicted branches one can take a semiclassical limit and find 
vacua corresponding to (classically) unbroken U(Ni) x U{N2) gauge groups. 

1 We assume that all Ni ^ 0. See Appendix C for some comments about the case when some 
of the Ni are zero. 



4 



1.2. Dessins From Gauge Theory 



Let us continue our analysis of the U(4) example. In Figure 2, we show a typical 
configuration of the zeroes of the polynomials that appear in (|1 . 2j ) , which correspond to 
points in the branches shown in Figure 1 near the semi-classical limits. The line segments 
that are drawn represent the branch cuts. Of course, it is not essential to draw the cuts 
through the zeroes of P(z) (denoted by o) but this will have a mathematical significance 
when we formally define a dessin. On the other hand, the cuts naturally pass through the 
zeroes of H$(z) (denoted by a bivalent • node in the drawing) since a small deformation 
away from the M = 1 branch will split the double zeroes. As is clear from the figure, these 
look like disconnected branchless trees. 





Figure 2: Zeroes of P&(z) (denoted by o), zeroes of Fa(z) (denoted by univalent •) and zeroes 
of H%(z) (bivalent •) near the semi-classical limits with Ni = N2 = 2 (left) and Ni = 1, N2 = 3 
(right). Edges represent branch cuts. 

It turns out that one of the two branches in Figure 1 has U(2) x £7(2) as the only 
semi-classical limit H. We now focus our attention on this branch. We further tune the 
parameters of the superpotential so that the corresponding M = 1 vacuum is an isolated 
singular point where Fl<i(z) develops a double root. There is only one such point in the 
branch we have chosen and it naturally leads to a connected graph. This is shown in Figure 
3. We have omitted the zeroes of Pa(z) so as to not clutter the figure. 

It turns out that such connected trees show up in the moduli space whenever the 
Seiberg-Witten curve develops isolated singularities. Examples of such "rigid curves" in- 
clude the maximally confining points |l|,|12| and the generalized Argyres-Douglas points 
fl3|1 . Such singularities arise when some of the zeroes of F 2n and Hj^_ n coincide, -F 2n (z) 
develops double roots, etc. 



U{2) x U(2) 




Figure 3: Evolution of the two branchless trees in one J\f = 1 branch of the U(4) gauge theory 
starting near the U(2) x U (2) semi-classical limit and reaching an isolated singularity where we 
get a connected tree. 

The connected trees that appear at such special points in the moduli space are precisely 
the "dessins d'enfants" that Grothendieck introduced as a tool to study the structure of 
the absolute Galois group Gal(Q/Q). This is the group of automorphisms of Q that leave 
Q fixed!. As mentioned earlier Gal(Q/Q) acts faithfully on the dessins. This means that 
the only element of the group that leaves every dessin invariant is the identity. The set 
of dessins is then partitioned into orbits under the action of the group. (We exhibit how 
Gal(Q/Q) acts on the dessins in Section 2.4.) One way of learning about the Galois group 
is to construct a complete set of invariants that distinguish dessins that belong to distinct 
Galois orbits. This is one of the goals in the study of dessins in mathematics and we will 
discuss some of the known Galois invariants in detail in section 3. 

In the following sections we show that the known order parameters that distinguish 
different branches of Af = 1 vacua can be thought of as Galois invariants. In particular we 
prove that the "confinement index" introduced in H is a Galois invariant and can be given 
a purely combinatorial interpretation. Interestingly, we will find that certain operations on 
the gauge theory side, such as the "multiplication map" introduced in |TtJ , have a precise 



interpretation as operations on the dessins. We believe that solving the non-rigid problem 
in ( |1.2j ) before specializing to singular points might lead to a new and useful perspective in 



2 Here Q is the field of rational numbers and Q is its algebraic closure. For more details see 
Appendix A. 
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the study of dessins d'enfants. Conversely we will also see that this correspondence leads 
to open questions regarding the interpretation of interesting mathematical invariants in 
the gauge theory. 

More explicitly, we would like to argue that the special points where the dessins make 
their appearance can be thought of as special phases of the Af = 1 gauge theory. This 
would imply the existence of appropriate order parameters which distinguish these special 
points from generic points in the M = 1 branch to which they belong. We provide some 
evidence for this in the examples in Section 5.2. 

So far we have seen in an example how dessins can arise at isolated singular points in 
the moduli space of a Seiberg-Witten curve. We will now show that given a dessin, one 
can associate to it a polynomial equation which corresponds to a singular Seiberg-Witten 
curve, of the type discussed in this section. It turns out that this is precisely equivalent to 
the content of the Grothendieck correspondence which we discuss next. 



2. The Grothendieck Correspondence 

2.1. Mathematical Preliminaries 

The Grothendieck correspondence is a bijection between classes of dessins and special 
classes of maps on punctured Riemann surfaces called Belyi maps. At first sight this might 
seem far removed from gauge theory physics but we will describe a precise route to arrive 
at the correspondence between Seiberg-Witten theory and the dessins d'enfants. 



The first ingredient in the correspondence is the Belyi map. A Belyi map |14[ is a 
holomorphic map from any punctured Riemann surface to P 1 with exactly three critical 
values at {0, l,oo}. In J7| Grothendieck showed that any dessin can be constructed from 
a Belyi map. For the purposes of our discussion, we will restrict to dessins drawn on a 
Riemann sphereR We show two simple examples in Figure 4. The key result that makes 
explicit the relation to Seiberg-Witten theory is that Belyi maps are obtained as solutions 
to certain polynomial equations. 



All the mathematical discussion in this section can be generalized to a general Riemann 



surface and we refer the interested reader to [15]. 





Figure 4: Examples of dessins; the bipartite structure of the graph is manifestly shown. 

As we will see, these equations have a natural interpretation as Seiberg-Witten curves 
that describe particular degenerations of Riemann surfaces, such as those we have already 
encountered in the introduction. More generally, we find that whenever the Seiberg-Witten 
curve factorizes so as to give rise to a "rigid curve" , i.e. without free parameters, one can 
associate a dessin to itR Our goal will be to set up a dictionary that maps the relevant 
quantities in mathematics, related to the Galois group action on dessins, into gauge theory 
language and vice versa. 

2.2. Dessins From Belyi Maps 

We now state without proof some basic mathematical facts which are crucial to estab- 
lish the relation between the dessins and Seiberg-Witten theory. See [|15| for a thorough 
discussion and for a full list of references. The main result we will use is the Grothendieck 
correspondence, which connects the theory of dessins with algebraic curves defined over 
Q, the algebraic closure of Q. Let us see how this comes about. 

Consider an algebraic curve X defined over C. Such a curve is defined over Q if and 
only if there exists a non-constant holomorphic function / : X — - > P 1 such that all its 
critical values lie in Q. A theorem by Belyi |14[ gives a very striking result: X is defined 



over Q if and only if there exists a holomorphic map / : X — > P 1 such that its critical 
values are {0, 1, oo}. 

A map [3 : X —>■ P 1 with all its critical values in {0, 1, oo} is therefore called a Belyi 
map. A Belyi map is called clean if all ramification degrees over 1 are exactly equal to 2. 
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This is up to a shift of z in (1.2). In gauge theory, this corresponds to the overall U (1) degree 



of freedom that decouples from the strong dynamics in the infrared. 
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Let us give a simple example that will be very relevant in the rest of the paper. Let 
X = P 1 and (3 a polynomial. To guarantee that all critical points that map to 1 have 
ramification degree 2, we set 

(3(z) = l-P 2 (z), (2.1) 

where P(z) is a polynomial. Let us see under which conditions f3 is a clean Belyi map. 
The critical points are computed as the zeroes of 

^- = -2P(z)P'(z). (2.2) 

This means that the zeroes of P(z) are critical points. Their ramification degree is 2 since 
P(z) is squared in f3 and their critical value, i.e., (3 evaluated at a zero of P(z), is 1. All 
we need is that the remaining critical points, which are precisely the roots of P'(z), have 
critical value 0. In other words, they must also be roots of 1 — P 2 (z). Up to the freedom to 
shift z, these conditions have only a discrete number of solutions. These are Belyi maps. 



o 
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Figure 5: We show how the dessin is the pie-image of the interval [0, 1] under the Belyi map. 
Note that as we move from to 1, the number of lines emanating from a given pre-image of is 
given by the ramification degree of the map at those points. Since the map is clean, exactly two 
lines meet at each pre-image of 1 . 

9 



We now have the ingredients to formally define a dessin: for the purposes of this 
paper we define a dessin d' enfant on the sphere as the pre-image under a clean Belyi map 
of the interval joining and 1 in P 1 . In other words, the dessin D associated to a clean 
Belyi map f3 is D = /3 _1 ([0, 1]) C X. We show this pictorially in Figure 5 for the case 
(3(z) = l-Pi(z)=F 4 (z)Hf(z). 

Such a dessin has a natural bipartite structure given by assigning a • to the preimages 
of and o to the preimages of 1. We will refer to a pre-image of as a vertex of the 
dessin. An edge is a line segment between two vertices that contain exactly one pre-image 
of 1. For example the second dessin in Figure 4 is clean (while the first is not) so that the 
notion of edges and vertices as we just defined makes sense: it has 7 edges and 7 vertices. 
In what follows, we will restrict our attention to clean dessins and refer to them simply as 
dessins. Likewise, the corresponding clean Belyi maps will be simply called Belyi maps. 

Also important for characterizing the dessins are the preimages of oo, denoted by x. 
There is one pre-image of oo for each open cell enclosed by a set of edges. For example, a 
dessin is a tree if and only if the preimage of oo is a single point. 

The study of dessins on the sphere is important because the absolute Galois group 
Gal(Q/Q) acts faithfully on them. As mentioned before, the absolute Galois group is 
the group of automorphisms of Q that leaves invariant Q and it is a remarkably complex 
object. We give a basic introduction to the Galois group in Appendix A. It can also be 
shown that not only is the action of Gal(Q/Q) on genus-0 dessins faithful, but so is the 
action on the much smaller set of trees. 

The main thing to take away from this section is that one can map the problem of 
classifying dessins to the problem of classifying Belyi maps f3. As we have seen, these are 
a special class of rational functions on the Riemann sphere that satisfy the conditions in 
the definition above. We now turn to the question of how Belyi maps corresponding to 
a given dessin can be explicitly constructed. This will naturally lead us to gauge theory 
physics. 

2.3. Belyi Maps From Polynomial Equations 

Consider a dessin D on the sphere with N edges. Let V = {ui, . . . , Uk} where Ui is 
the number of vertices (pre-images of 0) of valence i. We choose k to be the maximum 
vertex valence in D. Let C = {v\, . . . , v m } where V{ is the number of faces with % edges. 
Again we choose m to be the maximum face valence in D. The lists V and C are called 
the valency lists of D. 
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Let G Vi (z) and J Ui (z) be polynomials of degree Vi and Ui respectively, with undeter- 
mined coefficients. Take one polynomial for each element in V and in C. Pick i$ to be the 
valence whose vi is the smallest non zero value in C. Then let all polynomials G Vi (z) and 
J Ui (z) be monic except for G Vi (z) which we choose to be of the form 

G ViQ {z)=a(z"o- 1 + ...). (2.3) 

In other words, we have chosen the coefficient of z Vi o to vanish and we have factored out 
the coefficient of z Vi o _1 which we call a. 
Now construct the two polynomials 

k m 

M*) = n j ^ > b & = n ^ w ■ (2-4) 

j=i i=i 

Then if A(z) and B(z) are such that there exists a monic polynomial Pn(z), with iV equal 
to the number of edges of the dessin, and which satisfies the polynomial constraint 

A(z)-B(z)=P 2 N (z), (2.5) 

then 

M -1+*M-A& (26) 

P{Z) ~ + B{z) ~ B(z) [ } 

is a rational clean Belyi map. Note that the polynomial equation is rigid, in the sense that 

there are no coefficients in the polynomials which are free parameters!. There are thus 

only a finite number of solutions to ( |2.4j ). The discussion has been rather abstract so far, 

so let us illustrate the various concepts with some simple examples which have already 

been discussed in | fL6"f . 

Consider the dessin in Figure 6A which has 6 edges. From the figure, we see that 

each open cell is bounded by 3 edges and every vertex is trivalent. The valency lists are 

therefore of the form V = {0, 0, 4} and C = {0, 0, 4}. From the discussion above, we need 

A(z) = J%(z) and B(z) = G\{z) such that they satisfy the polynomial equation 

Pl{z)+Gl{z) = 4{z). (2.7) 



Up to a shift and rescaling of z. We will return to this point in Section 2.7. 
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Figure 6A: Dessin corresponding to the solution of the factorization problem ( |2, 7[ ), 

It turns out that there is only one solution to the polynomial equation ( [2.71) 01 ■ 
However, in general, such polynomial equations have more than one solution. For instance, 
the polynomial equation 

PUz) + Gl{z) = J!(z)Ji(z). (2.8) 

turns out to have two solutions |1(J. This is because for the same valency lists V = {0, 4, 4} 
and C = {0,0,0,0,4} (which we can infer from the polynomial equation), there are two 
dessins one can draw. These are shown below in Figure 6B. We shall revisit this specific 
example in Section 6 in much more detail. 




Figure 6B: Dessin corresponding to the solution of the factorization problem f |2.5|j . The bivalent 
• vertices are roots of J±{z) while the trivalent • nodes are the roots of Ja{z). 

The key result which we will use from now on is that for every dessin D with valency 
lists V and C there exists a solution to the factorization problem ( [2.50 such that the 
corresponding Belyi function gives D = /3 _1 ([0, 1]). This is a simplified version of the 
Grothendieck correspondence [0. 
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2.4- Action Of The Galois Group 

So far we have mentioned repeatedly that the absolute Galois group V = 
acts faithfully on dessins. We now show how T acts on the dessins via the Belyi map 
Let D be a dessin such that D = /3 _1 ([0, 1]). Furthermore, let f3 be of the form 

A(z) z 2N +a 1 z 2N - 1 + ... + a 2N 



B(z) z L + b 1 z L ~ 1 + ... + 6 7 



where A(z) and B(z) solve the polynomial equation (|2.5|), N is, as before, the degree of 
P(z) and L = YliLi ^ v i- Then V acts on D by twisting the coefhcientsil in (3 . For g £ T, 
D g is defined to be the dessin obtained by the action of g on D, i.e. D g = /3~ 1 ([0, 1]), 
where 

R = AM = z 2N +9(ai)z 2N - 1 + ... + g(a2N) . q) 

Pd B g (z) z I '+g(b 1 )zL-i + ... + g(b L ) ' K ' } 



Thus, given a solution to the polynomial problem (|2.5|) it is easy to understand how 
the dessin changes under the action of the Galois group. However, given two dessins, it is 
in general very difficult to tell whether they belong to the same Galois orbit or not. This is 
the central problem associated to the dessins d'enfants. Later we will discuss several Galois 
invariants that mathematicians have introduced in order to distinguish dessins that belong 
to distinct Galois orbits by studying the combinatorial data associated to each dessin. 

2.5. The Identification 

We have already discussed the Seiberg-Witten curves for pure gauge theories in ( |1.1| ) . 
Recall that for a U(N) gauge theory with L < 2N massive flavors with masses given by 
rrii, the Seiberg-Witten curve that captures the infrared dynamics of the gauge theory is 
the following hyperelliptic Riemann surface []l7|jr^1 : 



L 

y 2 = (det(zl - $)) 2 - 4A 2N ~ L Y[(z + m % ) . (2.10) 

i=i 

We would like to propose the following identification and argue that it is a useful 
one. Let us identify objects in (|2.5|) and in (|2.10|) as follows: Pn( z ) = (det(zl — <&)}, 



6 The a; (and bj) are algebraic numbers; i.e. they are solutions to some polynomial equations 
with coefficients in Q. The solutions to such equations include a, along with other algebraic 
numbers which are, by definition, in the Galois orbit of a%. Twisting by the relevant group element 
of r here refers to choosing another element in the orbit of a^. For a more formal discussion, refer 
to Appendix A. 
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B(z) = —4A 2N ~ L Y[f=i( z ~ m i)- I n particular, a = — 4A 2N ~ L . In each case, the precise 
form of A(z) in (|2.4j ) defines the special point in the Coulomb moduli space we are looking 
at. Since dessins are associated only to rigid factorizations of the Seiberg-Witten curves, 
they appear at isolated singular points in the moduli space of the Af — 2 gauge theory. 

At this point, it appears as if the relation to Seiberg-Witten curves is purely at a 
formal level. We will show in what follows that this is more than a superficial similarity 
and we exhibit features of the gauge theory that have a natural interpretation as operations 
on the dessin. For this we have to abandon our Af = 2 point of view and deform the theory 
to Af = 1 by a tree level superpotential as reviewed in Section 1.1. 

We mentioned earlier that the absolute Galois group acts faithfully on the set of all 
trees. For most part of the paper we will restrict our discussion to the set of trees and 
only in Section 6 will we discuss dessins with loops. 

2.6. Trees On The Riemann Sphere: Refined Valency Lists 

Since trees have only one open cell (with the associated vertex at infinity) the Seiberg- 
Witten curve associated to the dessin is that of the pure U (N) gauge theory: 

y 2 = P 2 N (z) - 4A 2N . (2.11) 

By itself, the curve in ( |2.11| ) does not correspond to any dessin, but if we tune the pa- 
rameters in Pn(z) so that we are at an isolated singularity in the moduli space, the curve 
factorizes, and the zeroes of the polynomials involved will describe vertices of a dessin. 
This also means that the Belyi map (|2.6|) is a polynomial 



P[Z) ~ B{z) ~ 4A 2N ' { } 

where Pn(z) solves the factorization 

k 

(P N (z) - 2A N )(P N (z) + 2A N ) = H(J Uj (z)y . (2.13) 

3=1 

From the expression it follows that the two factors on the left cannot have any factors in 
common. Thus, for trees, the problem always reduces to solving two lower order equations 
of the form 



P N (z) -2A N = l[(Q u -(z)y 



P N (z) +2A N = H(R u +(z)y 

3 = 1 
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such that it ■ + uf = it, for every j. One can now define a new bipartite structure on 
the dessin by assigning a +(— ) to every zero of R u + (Q u -) such that if a given vertex 

J 3 

(pre-image of 0) is of one sign, every one of its neighbors is of the opposite sign. The 
bipartite structure is unique up to an overall sign flip. This leads to a more refined valency 
list {V + , V~} than the {V, C} introduced earlier! where V± denotes the positive/negative 
valency list. For instance, from ( |2.14j) is the number of j-valent vertices of type "+" . 



With the refined valency list, the trees can be redrawn as shown in the figure below. 




Figure 7: Refined valency list for the trees. V + = {2, 0, 0, 1} and V = {3, 0, 1} for this case. 

Note that in the figure on the right we have removed the preimages of 1 depicted as 
o on the left. This is common practice in the literature when dealing with trees. There is 
one more common convention which is to depict elements in V + by • and elements in V~ 
by o; we have not adopted this convention in order to avoid confusion and we simply add 
± to the «'s as in Figure 7. 

Clearly, the dessins that arise from different valency lists belong to distinct Galois 
orbits; we will comment more about this in the next section. 

Interestingly enough, there is a related splitting of the polynomial equation in the 
gauge theory. In || while solving the non-rigid problem (|l]3|) it was found that the M = 1 
branches are classified by the integers where s± refers to the number of double 

roots in each of the factors on the left hand side of ( p.!4j) . This is already a hint that the 
mathematical goal of classifying dessins according to Galois orbits might be closely related 
to the more physical problem of studying the branches of M = 1 vacua in gauge theory. 
We will see this in more detail in Section 4. 



7 C contains just one element and is trivial for the case of trees. 
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2.7. Equivalence Classes Of Trees 



Let us consider the equations (|2.14|) in more detail. These equations have two free 
parameters corresponding to the scale and shift of z. In other words, if pjp(z) is a solution 
then pffi(z) = pjp(az + b) is also a solution. In physics as well as in mathematics, it is 
natural to consider monic polynomials. This means that a is restricted to be an N th root 
of unity. 



The trees constructed using the Belyi map (grp with P^ 1 > (z) and P ( N ' (z) are identical 



except for a displacement or rotation in the z plane. In the mathematical literature, such 
trees are considered equivalent and one considers equivalence classes of such Belyi maps. 
The Grothendieck correspondence is in fact an isomorphism between equivalence classes 
of Belyi maps! and children's drawings. 

Every tree has associated to it a number field (see Appendix A for a primer on field 
extensions) [J15| , determined by the field of definition of the polynomial Pn(z) that gives 



the corresponding Belyi map ( |2.12| ). This might be a little puzzling at first, since the 
transformation z — > az + b can in general involve arbitrary algebraic numbers. This means 
that the number field associated to trees that differ by translations and rotations can be 
different. On the other hand we have just said that such trees are taken to define an 
equivalence class on which Gal(Q/Q) acts. 

The resolution to this puzzle is that, although the Galois group acts nontrivially on all 
these trees, there is always a way of choosing the tree with the simplest number field ]19| 
as a representative of the equivalence class. It turns out that the action of Gal(Q/Q) on 
just the representatives of each class is faithful. Therefore, for the purposes of studying the 
absolute Galois group one uses the shift and scale of z to pick the simplest representative. 

All these statements have a counterpart in physics. The freedom to shift by b corre- 
sponds to the fact that the underlying theory is U(N), as opposed to SU(N). The overall 
£7(1) decouples in the IR and gives rise to this shift degree of freedom. Just like in mathe- 
matics, one can use this shift to bring any tree level superpotential to a form that displays 
the M = 1 branches, introduced in Section 1, most clearly. We will use this in section 5. 

More intriguing is the meaning of the rescaling by an TV-th root of unity. In physics, 
this corresponds to different kinds of confinement distinguished by the behavior of com- 
binations of Wilson and 't Hooft loop operators. Roughly speaking, the trivial root of 



The equivalence class of Belyi maps is, in fact, up to any SX(2,(D) transformation. However, 
we have used one of these to put the pole at oo. Thus only shift and scale transformations remain. 
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unity corresponds to usual confinement while the other roots correspond to oblique con- 
finement ||. Quite nicely, the mathematical criterion of choosing the simplest number 
field corresponds in physics to choosing the phase with usual confinement. 

It would be very interesting to explore the connection between the "not-so-simple 
number fields" and the oblique confining phases. However, since our goal is to establish 
a connection between dessins (in terms of equivalence classes) and gauge theory, we will 
restrict our study to physics phases with only usual confinement. 

2.8. Example: The Maximally Confining Af = 1 Vacua 

For now, let us discuss as an example the simplest tree one can draw: a branchless 
linear tree with N edges. Such a tree has 2 vertices of valence 1 and N — 1 vertices of 
valence 2. From the general discussion above, it is easy to write down the corresponding 
Seiberg-Witten curve for this case: 

P 2 N {z)-A={z 2 -4)H 2 N _ 1 (z). (2.15) 

Here we have set = 1. 

This Seiberg-Witten curve corresponds to points where N— 1 mutually local monopoles 
go massless. The M = 1 vacua are obtained by perturbing the M = 2 theory by a mass 
deformation, with Wtree = ^Tr <f> 2 . The condition = 1 has N different solutions that 
correspond to the N different maximally confining M = 1 vacua§. However, as discussed 
above, we take the simplest solution, i.e., A = 1. The solution to this equation is well 
known and given in terms of Chebyshev polynomials!^. In Figure 8 the zeroes of Pn(z) 
and Hn-i(z) have been depicted showing how the branchless tree arises. 

The importance of the branchless tree lies in the fact that these appear at the inter- 
section of the M = 1 branches, the point marked by a cross in Figure 1. They also appear 
near the semi-classical limits (A — > 0) as shown in Figure 2. Although these truncated 
branchless trees cannot be thought of as dessins, they seem to be "building blocks" that 
come together to create a dessin at an isolated singularity. We discuss this point in more 
detail in the conclusions. 

So far we have been rather loose in the language employed to discuss aspects of the 
factorization problems. Techniques from both the physics and mathematics literature have 

9 The reason for the name is that the low energy gauge group is just U(l) C U(N). 
10 Here the shift symmetry is used to set < Tr $ >= 0. 
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been used interchangeably. We now turn to a more systematic discussion of how the dessins 
are classified from a mathematical point of view. We will follow this up with a review of 
how the M = 1 vacua are classified from a physics point of view. 

O zeroes of Pn 

# zeroes of Hn-i 

# zeroes of z — 4 



N even/odd 



+ 



7+ +/- 



Figure 8: The dessin that corresponds to the maximally confining vacuum. Case by case, it is 



obtained by plotting the zeroes of the polynomials that solve the factorization problem ( 2.15 ). 
See for example Figure 17 for the N = 6 plot. 



3. Invariants 

3.1. Invariants From Mathematics 

One way of learning about the structure of the absolute Galois group Gal(Q/Q) is by 
constructing a complete set of invariants under the action of Gal(Q/Q) such that any two 
dessins that do not belong to the same Galois orbit will disagree in at least one invariant. 
Such a complete list of invariants is currently not known although many invariants have 
been constructed. In this section we review the most basic invariants0 associated to 
dessins that are related by the action of Gal(Q/Q) |15|, p0[ . 

• Valency Lists. The most intuitive invariants are the valency lists V and C introduced 
in Section 2. These are clearly invariants, since as we saw they are determined by the 
form of the polynomial equation which is invariant under the action of Gal(Q/Q) . It is 
sometimes possible to define more refined valency lists, corresponding to different ways of 
solving the polynomial equation. We have already seen this for the case of trees, where we 

11 More invariants than those discussed here are known, but for the purpose of our analysis we 
will concentrate on those that can be most easily computed explicitly. 
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introduced the {V + , V~} valency lists. As we discussed, this possibility of constructing a 
new invariant has a nice counterpart in physics; it will be further clarified in the examples 
that follow. Note that by concentrating on dessins coming from the same factorization 
problem we can forget about the valency list invariant since all dessins constructed this 
way have the same valence list. Therefore, the search is for other invariants that will 
distinguish different Galois orbits. 

•Monodromy Group. Every dessin is associated to a cover of P , defined by the Belyi 

map 

/3 : £ -> U = P 1 \ {0, 1, 00} , 

that maps the edges of the graph on the Riemann surface £ into the open segment 01 on 
P 1 . If the graph has iV edges, where we count the number of edges to be equal to the 
number of pre-images of 1, the map (3 gives a 2 N- fold cover of U (see Figure 9 for an 
illustration). 

Consider 7ti(U, 01) , the homotopy group of paths in U that begin and end on a point 
of 01 . Since a closed path based on 01 in P 1 can be mapped into a path between any two 
of the 2N segments in the fiber over 01 , any given element of iti(U, 01) acts on the dessin 
as a permutation of the half-edges (that go between a filled and unfilled vertex in Figure 
9). Therefore, the covering map (3 induces a map from iti(U, 01) to S2N ■ Let us denote by 
o"o the permutation corresponding to circling once the point z = on P 1 and by o~\ the 
permutation corresponding to circling once the point z = 1 . Recall that the dessins have 
a bipartite structure that keeps track of whether a vertex is mapped to or 1 by the Belyi 
map. One can convince oneself that oq is the element of S2N that permutes cyclically the 
edges incident on each vertex (that maps to 0) and, similarly, <j\ is a cyclic permutation 
of the edges incident on each pre-image of 1. The subgroup of S2N generated by o"o and 



o"i is the monodromy group of the dessin and it is a Galois invariant ^21 . 

Let us consider the example of the dessin in Figure 9. There are 8 half-edges. The 
monodromy group is generated by the following permutations: 

a = (1,7, 6)(2,3)(4,5), 

(3.1) 

(7! = (1,2) (3, 4) (5, 6) (7, 8). 
We will give the explicit monodromy groups for the examples we will encounter later. 
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• o 

1 

Figure 9: Left: Example of a dessin with N = 4 edges. Right: The 8 half-edges come from the 
preimage of the open interval 01, i.e. the preimage of and 1 are not shown. 

•Belyi Extending Maps. It is possible to construct new invariants by composing the 
Belyi map of interest with any Belyi-extending map and then computing the valency lists 
or the monodromy group of the new dessins obtained this way [p0 |. A Belyi extending 



map is a Belyi map a : P 1 -> P 1 defined over Q , such that its composition with any other 
Belyi map does not change the associated number field. For our purposes, we relax the 
definition slightly: by a Belyi extending map here we will mean any map a defined over Q 
that can be composed with a Belyi map f3 to give another Belyi map (3 a := a o (3. 

Let / be any invariant of a dessin D a , where D a = /3~ 1 ([0, 1]); the claim |20| is that 
/ is also an invariant of the dessin D = /9 _1 ([0, 1]). For example, for the Belyi extending 
map ot2{z) = Az(l — z) , the monodromy group of D a is the cartographic group of D, which 
is known to be another Galois invariant. Later we will prove that the multiplication map 



of []lQ|Jq1 is the physical realization of the Belyi extending map «2- 



3.2. Invariants From Physics 

As discussed in the introduction, one problem that is very similar to the classification 
of dessins using Galois invariants is the problem of classifying branches of M = 1 vacua 
using order parameters, such as Wilson and 't Hooft loops. We will consider those N = 1 
vacua that are obtained in the infrared by starting with an M = 2 gauge theory and 
adding a tree level superpotential Wt Tee for the adjoint scalar The discussion of the 
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gauge theory order parameters in this section will closely follow that of [Q. In fact, what 
follows is just a summary. We refer the reader to || for all relevant details and proofs. 

• Confinement Index. In a U(N) gauge theory a natural order parameter is the 
expectation value of a Wilson loop W in, say, the fundamental representation. The Wilson 
loop in the tensor product of r fundamental representations is W r . Clearly, for r = N there 
is no area law, for it is equivalent to a singlet representation (due to electric screening). 
A measure of confinement is the smallest value of r, which can only be between 1 and 
N, for which W r does not exhibit an area law. Such a value is denoted by t and it is 
called the confinement index. When the gauge group is broken classically to a product of 
factors U(Ni) x U{N2) X ... x U(N n ) one has to also use the 't Hooft loop H to determine 
the confinement index. By embedding a 't Hooft-Polyakov magnetic monopole of the full 
U(N) theory in any two of the U(Ni) factors, we take into account magnetic screening. 
If for each U(Ni) we get that W[ i Hi has no area law, this implies that in the full U(N) 
theory W Ti ~ rj has no area law. The relative sign comes from the fact that the magnetic 
monopole sits in both groups with opposite charges. 

Therefore, after taking into account electric and magnetic screening, the confinement 
index is given by the greatest common divisor of the Ni and bi = ri — r^i. These two sets 
of quantities, N^s and fej's, will have a very clear combinatorial meaning, which will allow 
us to compute the confinement index just by inspection of any dessin. 

As a preparation for that let us mention that both set of quantities are encoded in 
the expectation values of the generating function for chiral operators Tr $ s , given by 



It turns out that the periods of T(z)dz, thought of as a meromorphic differential on y 2 = 
W{ iee (z) 2 + f n -i(z), are related to the iVj's and b^s as follows: the iVj's are the periods 
of T(z)dz on the A-cycles and the b^s are the periods of T(z)dz on the S-cycles (for 
an appropriate choice of basis). The b^s measure the relative theta angle of U(Ni) and 
U(Ni+i). Moreover, one can show that in the M = 1 branch with confinement index 
t, T(z)dz = tT(z)dz, where T(z) dz is the generating function for chiral operators in 
a Coulomb vacuum (which have t = 1) of a U(+f) theory ||. The fact that the two 
generating functions are related by a multiplication by t has an important consequence: 
all confining vacua with confinement index t are obtained from Coulomb vacua by using 



the multiplication map by t [10. qj. The definition and discussion of the multiplication map 
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is given in Appendix B. We also discuss this further in Section 4.1 where, for the specific 
case of t = 2, the multiplication map will be shown to coincide with the Belyi extending 
map «2 of [15]. 

• Holomorphic Invariants. In cases when the rank of the low energy gauge group is 
too high, i.e. when the degree of the tree level superpotential is large, there is always at 
least one iVj which is equal to 1. The precise condition is degW(z) > N/2: when this 
condition is satisfied, the confinement index is always one. One might naively think that 
there is only one branch since we have exhausted the standard order parameters. However, 
it is possible to show that there are many branches, all of them having Coulomb vacua. 
This is the problem that motivated the search for non-conventional order parameters in 
||. It turns out that the discussion that follows also applies for superpotentials of any 
degree. 

The new non-conventional order parameters proposed in || are obtained by studying 
relations between the vacuum expectation values of different chiral operators that can 
be defined in the theory. The expectation values of chiral operators become holomophic 
functions on the moduli space of vacua due to supersymmetry. It turns out that, at least 
in the examples considered in H, these functions satisfy different polynomial constraints 
in different branches. The problem of whether the existence of these relations was the 
reason for the existence of the different branches or viceversa was left as an open question. 
For the purposes of this paper, we take the former as the correct point of view. In fact, 
we will see that in the mathematical literature the refined valency list for trees gives very 
similar information as the relations between holomorphic functions found in ||. 

The chiral operators of relevance are Tr <& r W a W a , whose appropriately normalized 
vacuum expectation value is denoted by t r = — (l/327r 2 )(Tr & r W a W a ) . In terms of the 
reduced Seiberg-Witten curve y 2 = F 2n (z), they can be computed as0 



The relations introduced in |J to distinguish between different branches are polynomial 
equations in the t r 's. 

12 Strictly speaking, the curve needed for the computation is given by the matrix model curve 
of the Dijkgraaf-Vafa correspondence. However, in cases when none of the Ni are zero, the matrix 
model curve is the same as the reduced Seiberg-Witten curve. 




(3.3) 
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The different branches that these relations distinguish are determined by the dis- 
tribution of double zeroes of the curve (2.6) in the two factors of P^(z) — 4A 2Ar , 
i.e. by the pair (s+,s_). In other words, if we start with the factorization problem 
Pl{z) - AK 2N = F(z)H 2 (z), we get 

P N (z) - 2A N = R N _ 2s _(z)H 2 s (z), 

(3.4) 

P N {z) + 2A N = R N _ 2s+ (z)H 2 + (z). 
In order to derive the relations it is convenient to write 



f \ fa Tfa 7~\ H a+ (z)R N - 2a+ (z) 4A" 

y{ Z ) = ^ R N _ 2s _{z)R N _ 2s+ {z) = ~— — Wl - 2 , , T-r . (3.5) 

v H s _(z) y J1 s+ {z)Rn-2s + {z) 

Since the integral defining the t r 's is around infinity, the computation can be carried out by 
expanding the square root. It is easy to see that if < r < + — 2 then only the leading 
term in the expansion contributes. It turns out that in order to distinguish different values 
of (s_|_, S-) all that is needed are relations between those (restricted) t r 's. Using the fact 
that A does not appear, by matching dimensions (which for t r is 3 + r) and by matching the 
charge under the £7(1)$ symmetry (which for t r is r), one concludes that the polynomials 
must be homogeneous in the number of $'s and the number of W a W a . Consider for 
example the case when s_ = 1 and s + = 3. Then one can show that tot 2 — 1\ = 0. 

We will see in the next section that (s+, s_) gives some information about the refined 
valency list of a dessin. However, the refined valency list contains more information. In 
section 5.2 we will show by means of examples that the extra information of the refined 
valency list can be obtained if one keeps the next to leading order term in the expansion 
of the square root. In other words, in the expansion 



= 1 - , # — + 0(A 2N /z 2N ) (3.6) 



H 2 + (z)R N . 2s+ (z) H 2 + (z)R N . 2s+ (z) 

the first term gives information about (s+, s_) while the second encodes the whole refined 
valency list. It would be interesting to find a combinatorial meaning of the higher order 
terms. It is important to mention that the extra relations we find distinguish the isolated 
point where the dessin appears from its neighbouring points in the Af = 1 branch. 
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4. Cross Fertilization 

Before we discuss some examples to illustrate our ideas, we would like to exhibit part 
of the dictionary between the mathematical and physical descriptions. First we show how 
the multiplication map [Til can be interpreted as an example of a Belyi-extending map 



2CH ; we also show that the information about the refined valency list of trees, described in 
Section 2.7, can be recovered by studying the holomorphic invariants. Most importantly, 
we give a combinatorial interpretation of the confinement index introduced in || . We then 
speculate on the relation between the classification of dessins and the study of phases of 
gauge theories in four dimensions and formulate a few precise conjectures. 

4-1. Multiplication Map As A Belyi Extending Map 

Consider a tree with N edges. The Belyi map has the form ( f2.12| ) 

The dessin corresponds to a rigid polynomial equation, which is a special point in the 
parameter space of the non-rigid factorization problem 

P 2 N {z)-AK 2N = F 2n (z)H 2 N _ n (z). (4.1) 

The multiplication map |Hj (with multiplication factor 2) guarantees that, if Pn(z) satis- 
fies the non-rigid factorization equation ( |4.1[ ), one solution to the factorization equation 

Pi N (z) - 4A 4N = F 2n (z)H 2 2N _ n (z) 



is given by (see Appendix B for details) 



P N (z) 



P 2N (z)=2A 2N T 2 i [ ^j , (4.2) 

where T 2 (x) — 2x 2 — 1 is a Chebyshev polynomial of the first kind. Since Pn(z) gives rise 
to a Belyi map, P 2 n(z) as defined in (|4.2| ) also leads to a Belyi map (3 2 n whose inverse 
image of the interval [0, 1] leads to a dessin with 2N edges. Therefore, by applying the 
multiplication map, we get a new Belyi map of the form 

8 -1 P ^-1 (0^ 1^ 
~ 1 - ~4A4N- -1-^2 -^n - 1 

_P%W L PM*)\ _ 4B (1 a . (4-3) 



a 2 o 13 



N 
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with 0:2(2/) = 4y(l — 2/)- We thus find that this map coincides with the Belyi extending 
map mentioned in Section 3.1, which relates the monodromy group to the cartographical 
group. 



In [|20| the author gives a prescription to draw the dessin associated with any Belyi 
extending map starting from the original dessin. Roughly, the procedure consists in draw- 
ing on P 1 the preimage through the Belyi extending map of the [0, 1] segment; then, one 
substitutes this new drawing in place of each segment of the original dessin. Let us apply 
this to our example. The map 02 is of degree two, so it covers [0, 1] twice and we expect it 
to double the number of edges of the dessin. More precisely, the critical point y = 1/2 is 
mapped by 02 to the vertex z = 1 of the [0, 1] segment. Therefore, we infer the following 
rule to draw the dessin obtained through the Belyi extending map 02 : 



One could check that for generic t the multiplication map by t is a Belyi extending 
map that substitutes to each edge in the original graph a branchless tree of length t . 
Naively, the dessins with confinement index 2 or higher would appear to be "scaled up" 
versions of smaller dessins. Indeed this is what one gets if one applies the multiplication 
map to the rigid factorizations that lead to the dessins. 

However, from the gauge theory point of view, one can also apply the multiplication 
map to the non-rigid problem ( |4.1| ) and then impose the constraints that leads to a rigid 
factorization problem. In other words, if J-'r/nr is the set of rigid/non-rigid factorizations 
and if M is the multiplication map acting on the factorizations J 7 , 



M(F R )C ^M(F NR )j , (4.4) 

where the last subscript R indicates that the factorizations are restricted to be rigid. This 
shows that the multiplication map is more than an operation to get new Belyi maps from 
old ones. Starting from a point in the moduli space of a U(N) gauge theory which is not an 
isolated singularity (so that there is no associated dessin), one can apply the multiplication 
map by t and sometimes obtain a singular point in the moduli space of the U(tN) theory 
where one can obtain a dessin. We will discuss such examples in Section 5. Moreoever, we 
will also see in Section 4.3, that applying the multiplication map to non-rigid factorizations 
is what allows us to prove that the confinement index is a Galois invariant. 
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4-2. Refined Valency Lists From (Refined) Holomorphic Invariants 



We have already mentioned in Section 2.7 that while solving the non-rigid problem 
( |1.2| ) it is convenient to classify the solutions in terms of integers (s+,s_), where s± refers 
to the number of double roots in either of the two factors (P(z) ± 2A N ). We would now 
like to relate these numbers to the refined valency list introduced in Section 2.7. From the 
definitions one can check that a given dessin with refined valency lists V + = {u~£} and 
V~ = } will appear in an J\f = 1 branch whose (s+, s_) values are given byi 

oo 

s± = J2 ku ik- ( 4 - 5 ) 

k=i 

Recall that Uf~ is the number of /c-valent vertices in the dessin. Clearly, the refined valency 
list contains more information than just the values of s±. 

From the discussion in Section 3.2, we have seen that the set of relations between the 
expectation values of chiral operators t r 's depends upon the distribution of double roots 
s±, i.e. different branches are defined by the different polynomial relations between the 
£ r 's. These chiral ring relations are satisfied at any generic point on that branch. However 
the dessins appear only at special points in that moduli space. A simple counting of 
parameters shows that at such points there will be more relations that are not generically 
satisfied. We will, in the examples to be discussed in Section 5, write down these extra 
relations satisfied by the t r 's explicitly. 

Since the generic relations seem to distinguish branches of M = 1 vacua, it is tempting 
to conjecture that these special points are isolated phases. In other words, they are distinct 
phases from their neighbours in the A/" = 1 branch. That this is the case is easy to see 
in some cases where the corresponding Af = 1 theory becomes superconformal in the 
IR. Moreover it is believed that there is always a choice of superpotential for which the 
resulting M = 1 vacuum flows to an interacting superconformal theory [|S|,|23,2~3 . 



13 In assigning a refined valency list to a dessin, there is an overall choice of sign in assigning 
+/— to the vertices. It follows that this amounts to exchanging s+ and S-. The same choice is 
also present in gauge theory. 
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4-3. Confinement Index As A Galois Invariant 

The physical interpretation of the confinement index t was given earlier in this section. 
We also explained how t can be computed from the periods of a particular meromorphic 
differential T(z)dz on the Seiberg-Witten curve. Here we first give a purely combinatorial 
description of the confinement index t. Then, we show that this is indeed a Galois invariant. 

Consider a given tree T, constructed as T= /3 _1 ([0, 1]) under a clean Belyi map (3(z) = 
1 — P(z) 2 , where P(z) is a polynomial. Let us concentrate on the preimages of under (3. 
There can be vertices with any valency. In particular, there must be vertices with valence 
one; this is a simple consequence of the fact that the Belyi map is clean. 

The procedure for computing t is the following: circle all vertices with odd valence. 
Choose any of the circled univalent vertices as the starting point. Move from the chosen 
vertex to the next circled vertex, say going clockwise around the tree, and count the number 
of edges between the two circled vertices; call it h\. Move from the second circled vertex 
to the next circled vertex. Again count the number of edges between the two vertices; call 
it hi- The most important rule to apply when going around the tree is that each circled 
vertex can be used only once as a starting point and only once as an end point. Therefore, 
if the next vertex was used previously both as an end point and as a starting point, one 
should skip it and go to the next one. Continue around the tree until there are no more 
unused circled vertices. The prescription makes sense because there is an even number of 
odd-valent vertices^. 

After completing this procedure one is left with a list of integers C = {hi, h 2 , . . . , hf}, 
where / is the number of odd vertices in the dessin. Then t is simply given by the greatest 
common divisor of the elements of C. Two simple examples are shown below. 




14 That there is an even number of odd vertices is clear from the fact that f3 is a polynomial of 
even degree. 
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This rule works because the vertices with odd valence are precisely the points between 
which one would draw the cuts in the gauge theory approach. Then, this definition can 
be seen to coincide with the definition of the confinement index introduced in ||. This 
follows from the fact that the integral of T(z) dz between successive vertices is 1. 

Having given a purely combinatorial definition of t we proceed to show that this 
is indeed a Galois invariant. The proof involves concepts and terminology reviewed in 
Appendix A. 

Let us consider all possible trees with a fixed number of edges N. If N is prime then 
the only possible values of t are t = 1 and t = N. The only tree with t = N is the 
branchless tree which is a tree defined over Q and hence it is its own Galois orbitM All 
other trees have t = 1 , so there is nothing further to prove. 

Consider an N which is not prime. Let N = p^ 1 . . . p r s s be the prime decomposition 
of N. Take any pi and consider the auxiliary polynomial P Pi (z) . Use the multiplication 
map by m = N/pi to produce what we called a non-rigid curve (for more details on the 
multiplication map see appendix B) 

1 -Tl (P w (z)) = (z))t£_i(J* (*))• (4-6) 

This depends on the Pi + 1 coefficients of P Pi (z). As discussed in section 2.2, the new 
Pn(z) = T m (P p .(z)) gives rise to a Belyi map 0{z) = 1 — P^(z) if and only if P' N (z) 
divides the right hand side of Q4.6Q . This is equivalent to imposing that the right hand 
side of ( |4.(j| ) has only N — 1 distinct roots. These conditions will give rise to polynomial 
equations for the coefficients of P Pi (z) = aoz Pi + . . . + a Pi+ i. Let the set of polynomials that 
must vanish be S = {/i(a), . . . , fj(a)}. Since all Chebyshev polynomials T m (z), U m -i(z) 
have coefficients in Q it follows that f(a) G Q[ao, . . . , a Pi+ i]. Therefore, there is a splitting 
number field Ks associated to the set S. It is a finite normal extension of Q and hence is 
left invariant by Gal(Q/Q). This means that the solutions form full Galois orbits. From 
the relation between the multiplication map and the confinement index t it follows that all 
such orbits can only have values of t which are multiples of m that divide N. This means 
that either t = m = N/pi or t = pim = N. As mentioned above, there is a single tree with 
t = N and therefore all other orbits must have the same value t = N/pi. 

15 Here what we have in mind is the branchless tree with the simplest number field, which in 
this case it is Q. 
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Consider now k = piPj and the polynomial Pk(z). As before, use the multiplication 
map by m = N/k. Following the same procedure we conclude that dessins arising this way 
can only have values of t which are multiples of m and that divide N. The only possibilities 
are t = N/k, N/pi, N/pj, N. As before, t = N gives a single tree. We have already proven 
that dessins with t = N/pi or t = N/pj can only come in full Galois orbits. Therefore the 
remaining dessins with t = N/k also arise in full Galois orbits. 

One can continue this argument by induction and conclude that any two dessins in 
the same Galois orbit must have the same confinement index t. Thus, the confinement 
index is a Galois invariant. 

4-4- Speculations About Dessins And Gauge Theory: Weak And Strong Conjectures 

We now have all the ingredients we need to formulate our conjectures precisely. We 
have already seen that there exist Af = 1 branches in pure gauge theory that are classified 
by order parameters such as the confinement index. Also, as mentioned in Section 4.2 
(and as we will show in some simple examples), at the special points where the dessins 
appear one has extra chiral ring relations. For an appropriately chosen superpotential, 
these points are believed to give rise to superconformal Af = 1 theories in the IR and 
thus define new phases. However, it is known that for some superpotentials, these theories 
might not be singular. This does not exclude the possibility that these points might be 
new phases, perhaps distinguished by less exotic behavior, such as extra massless states or 
smaller rank of the gauge group. 

Given the earlier discussion regarding the theory of dessins and the phases of super- 
symmetric gauge theory, our first conjecture should be fairly well motivated: all points 
where dessins appear correspond to special phases embedded in the Af = 1 branches which 
we call "isolated phases" . 

Our second conjecture, relating the phases to Galois orbits of dessins has a weak and 
a strong form. The strong form is easily stated: all Galois invariants are physical order 
parameters that can be used to distinguish the isolated phases of supersymmetric gauge 
theory with a given superpotential. 

On the other hand, based on the examples we work out in Section 5, we find that all 
the order parameters used to distinguish branches of gauge theories that meet a particular 
U(l) branch, defined below, are Galois invariants. This is what we refer to as the weak 
form of the conjecture. 
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Some comments are in order. By a 17(1) branch we mean a branch where the low 
energy group is a single U(l) C U(N). This can be called the maximally confining branch. 
This branch has the same dimension as the other branches JO. Any generic branch meets 
these U(l) branches at points where the corresponding dessin is a branchless tree with N 
edges. 

We believe that the weak form of the second conjecture is very likely to be correct and 
we provide evidence for it in the examples. The strong form is on much less firm ground. 
In particular, it relies on the correctness of the first conjecture and on assumptions that 
require much more study. 

It would be very important to gather more evidence for the strong form since, if 
true, it provides a striking connection between Grothendieck's program of unveiling the 
structure of Gal(Q/Q) via its action on dessins and the physics problem of classifying 
phases of supersymmetric gauge theories. Some of the most striking consequences would 
be for gauge theories with matter where physics order parameters are scarce. Almost all 
known Galois invariants would become new gauge theory order parameters. In the next 
sections, we will provide some evidence in support of these conjectures. 

A More General M = 1 Viewpoint And A Global M = 2 Viewpoint 

The possibility that there is always an "extremal" superpotential for which the M = 1 
U(N) gauge theory at one of the isolated singular points becomes superconformal in the 
IR motivates the following point of view. Up to now we have studied theories with a 
superpotential of a given degree. However, if we fix U(N) and vary the degree of the 
superpotential, the theories arising at the points where a dessin D appears can go from 
being non singular to singular in the IR. Let us denote by d(D) the smallest degree of an 
extremal superpotential for D. It is tempting to conjecture that two theories that arise 
at points corresponding to two dessins D and D' in the same Galois orbit will necessarily 
have d(D) = d(D'). In other words, d(D) might be a Galois invariant. We leave this 
problem as an interesting direction for future work. 

Finally, let us comment on yet another point of view. Suppose that we set the tree 
level superpotential to zero. Then we recover an M = 2 gauge theory. The Seiberg-Witten 
curve of section 2.5 becomes the curve describing the physics in the moduli space of vacua 
of a single theory. The valency lists C and V of section 2.3 simply encode information 
about the masses of particles in the theory. More explicitly, C determines the distribution 
of masses of fundamental hypermultiplets. V determines, up to modular transformations, 
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the charges of the various monopoles and dyons that are massless in addition to the U(1) N 
vector multiplets present at generic points. According to the cases studied in the literature, 
there is reasonable evidence to suspect that all such points are Af = 2 superconformal field 
theories^. 

It would be very interesting to explore the relation between the classification of dessins 
into Galois orbits and the classification of such Af = 2 superconformal field theories. A 
natural possibility, worth exploring, is that field theories giving rise to dessins in the same 
Galois orbit might be dual theories in some sense. 



5. Examples: U(6) Pure Gauge Theory 

In this section we would like to illustrate by means of examples the concepts we have 
covered up to now. We consider pure Af = 2 U(6) gauge theory broken to Af = 1 by a tree 
level superpotential. From the general discussion about dessins and polynomial equations, 
it follows that the Af = 2 moduli space contains an isolated singularity for every connected 
tree with 6 edges. It turns out that all such dessins can be obtained by just using a quartic 
superpotential. We discuss why higher degree superpotentials are not needed and list all 
dessins with their corresponding factorization problems in Appendix C. 

Here, we restrict our study to a cubic superpotential. The dessins we obtain are, of 
course, a subset of the general quartic superpotential but they turn out to exhibit all the 
relevant points of the physics-mathematics dictionary we have established. All U(N) gauge 
theories with N = 2, . . . , 6 were studied in detail in where one parameter solutions to 
the factorization problem 

P N {zf - 4A 27V = F A {z) H 2 N _ 2 (z) (5.1) 

are listed. In fact the analysis in this section can be easily repeated for all these cases. 
We choose U (6) because it is the simplest case that exhibits four different values of the 
confinement index, i.e. t = 1,2, 3, 6. 

The rigid factorizations corresponding to dessins are obtained by imposing suitable 
conditions on the solutions found in H along the lines we described in the introduction. 

16 Of course, the point with N — 1 mutually local massless monopoles is not a superconformal 
field theory. In this case one can write down, using S— duality, a local lagrangian describing the 
full behavior of the theory in the IR. 



31 



In Sections 5.1 and 5.2 we classify the dessins according to the M = 1 branches to which 
they belong and specify the order parameters that distinguish the special points where the 
dessins appear as isolated phases. 

From a mathematical point of view, in order to find the explicit Belyi maps, it is not 



necessary to start from the non-rigid problem (|5.1| ) . Instead, one solves the rigid problems 
directly. In Section 5.3 we will present the solution to all possible rigid factorization 
problems that can be derived from (|5.1|) along the lines of Jl6j by using differentiation 
tricks. This analysis shows explicitly the classification of trees into distinct Galois orbits. 
In Section 5.4 we will reproduce the same classification of dessins using some of the known 
Galois invariants. We will find that this parallels the classification of phases in gauge 
theory. 

We mentioned in Section 2.7 that, both in the physics and mathematical analysis, 
there is the freedom to shift and scale the z variable. From a physics perspective, it is 
natural || to shift the z variable in order to bring the superpotential to the canonical form 
W( iee (z) = z 2 — A and then analyze the M = 1 branches obtained by varying A. However, 
since our primary goal is to exhibit the dessins and where they appear in the gauge theory 
moduli space, in the examples that follow, we have followed the mathematical strategy to 
shift and scale z to put the solution in the simplest form possible, so that the number field 
associated to the tree is the simplest. 

5.1. U(6) Gauge Theory: A Physicist's Point Of View 

Let us now review the solution of || in detail. We then specialize to rigid factorizations 
by tuning the one free parameter available in the solutions to ( |5 . 1| ) . 

As described in H one can, first of all, classify M = 1 branches by the number of 
double roots in either factor (Pq(z) ± 2A 6 ). If we denote the number of double roots in 
either factor as (s+, s_), in our case, these can take the values (3, 1), (1, 3) and (2, 2). All 
these branches meet at vacua that have (s+, s_) = (3, 2) or (2, 3) at which the branchless 
tree discussed in Section 2.7 appears. The M = 1 branches are further classified by the 
confinement index and the non-conventional order parameters which we defined earlier in 
Section 3.2. Let us consider each value of (s+, s_) in turn. 

• The (3,1) And (1,3) Confining Vacua 

The general factorization problem (54.) is solved by the polynomials^ 



17 Refer to equation (3.47) in 
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P 6 {z) + 277A 6 = ((z - a) 2 (z -b)- 2eA 3 ) = R 2 (z) 

P 6 {z) -2 V A 6 = (z-a) 2 (z-b)((z-a) 2 (z-b)-4eA s ) = S(z) ^'^ 

with i] 2 = 1 and e 2 = r\. These polynomials can by obtained by the "multiplication by 2" 
map acting on either of the polynomials Pz[z) = [z — a) 2 (z — b) =F 2Ajj: 

PaW = 2AVT 2 (|W), (5.3) 

with Aq = k 2 A 6 , k 4 = 1 and e = ±k. The various signs and phases in these expressions 
are crucial so that all the M = 1 vacua are taken into account. However, as discussed in 
Section 2.7, since these lead to trees in the same equivalence class, we will drop such phase 
factors in what follows. 

The Rigid Quartic Factorization 

If we require that the polynomial R(z) has a double root (which is to say that its 
discriminant vanishes) we get the rigid factorization 

P 6 (z) 2 - 4A 12 = F A (z) H 2 {z) Qi(z) . (5.4) 

This fixes A to be 

A 3 = ^(a-b) 3 . (5.5) 



Substituting this into the polynomials in (|5.2|) and using the shift and scaling symmetry 
to set a = — 5 and b = 10, we get 



R\ Z ) = ( z _ 5 )4(^ + 10 )2 5 

S(z) = (z + 5) 2 (z - 10) (z 3 - 75z + 750) 



(5.6) 



From this, we see that the polynomials that solve the equation 

Pi(z) - 4(250) 4 = F 4 (z) H 2 (z) Qf(z) (5.7) 

are given by 



Q x (z) = z-h, H 2 (z) = (z + 5)(z + 10), F 4 = (z - 10)(z 3 - 75z + 750) , 
p 6 ( z ) = z 6 - I50z 4 + 500z 3 + 5625z 2 - 37500z - 62500 . 



(5.* 
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Plotting the zeroes of the polynomials leads to the tree in Figure 10. Let us make a few 
comments about the solution. All polynomials are defined over Q. From the discussion in 
Section 2.4 about the action of the Galois group, we see that the tree is left invariant; in 
other words, it is the only element in its Galois orbit. From (|5.3| ) we see that the tree has 
confinement index 2. However, observe that the tree is not a scaled up version of a smaller 
tree. This illustrates the point made in Section 4.1 and especially equation ( fOj) . One can 
moreover check that the combinatorial method for computing the confinement index, as 
explained in Section 4.2, also gives the correct answer t = 2. 

P 6 (2) 2 -4A 12 =F 4 (z) J ff|(z)QtW 

• roots of F4 
roots of Ho 

Z 6 - 

• root of Qi 
O roots of P 6 



-10 -5 1 10 



-6 



Figure 10: The tree obtained in the (3, 1) confining branch with a four-valent vertex. It has 
confinement index t = 2. 

The Rigid Cubic Factorization 

Starting from (|5.1|) one can also get another rigid factorization by tuning one of the 
zeroes of F±(z) to coincide with a zero of H^(z) : 

P*(z)-4A 12 =F 3 (z)Hl(z)Ql(z). (5.9) 

This leads to the condition a = b in (|5.2|) . This implies that S(z) has a cubic root at z = a. 
Using the shift and scale symmetry, we can set a = and A = 1. This leads to a very 
simple solution of ( |5T9| ) : 



P 6 ( z ) = z 6 ~4z 3 + 2, H 3 (z) = z 3 -2, F 3 (z)=z s -4 and Q x {z) = z . 
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The tree associated to the factorization is drawn below in Figure 11. 



Im z 




P*(z)-4A 12 = F 3 (z)H 2 3 (z)Ql(z) 

# roots of F3 

# roots of H 3 

# root of Qi 
O roots of Pq 



-o-» 

1.5 



Re z 



0.5 



Figure 11: Tree obtained in the (3, 1) confining branch with a trivalent vertex. It is a scaled up 
dessin, obtained by applying the multiplication map on a smaller dessin with 3 edges. 

Note that, unlike the quartic case, the multiplication by 2 is easily understood: this 
particular solution can also be obtained by first solving the rigid factorization problem 



P*(z)-4A e = Ql(z)F 3 (z) 



(5.10) 



and then applying to the resulting solution the multiplication map. The solution is once 
again defined over Q and the tree in Figure 11 is the lone element in its Galois orbit. 
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• The (2,2) Confining Vacua 

In this sector, the factorization problem (5.4) is solved by the polynomials^ 



P 6 (z) + 2A 6 = (z 2 + g - A 2 ) 2 (z 2 + g + 2A 2 ) , 

(5.11) 

P 6 (z) - 2A 6 = (z 2 + g + A 2 ) 2 (z 2 + g- 2A 2 ) . 

These polynomials are obtained by the "multiplication by 3" map; modulo phase factors, 
P&(z) in (|5.11| ) is given in terms of Pziz) = z 2 + g as 

P 6 (z)=2A G T 3 (^-) . (5.12) 



The trees in this branch will therefore have confinement index 3 
The Rigid Quartic Factorization 



One can set g = A 2 to get the quartic factorization flOl ), while one can rescale z to 
set A = 1. The resulting polynomials that solve ( |5.4|) are 

P 6 (z) =z A (z 2 + 3)-2, H 2 (z) = (z 2 + 2) , 

(5.13) 

Q 1 (z)=z and F 4 (z) = (z 2 + 3)(z 2 - 1) . 

Plotting the roots of the polynomials leads, this time, to the tree in Figure 12. The 
polynomials are defined over Q and so the tree is the only element in its Galois orbit. 

The Rigid Cubic Factorization 

Note that it is not possible to get the cubic factorization equation ( |5.9| ) by tuning 
the available free parameter. Thus, we do not find any trivalent tree in this branch of the 
moduli space. One can also show this using the combinatorial definition of the confinement 
index by trying to construct a trivalent tree with six edges and t = 3. 



18 Refer equation (3.50) in @. 
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P 6 (z) 2 -AA 12 = F 4 (z)Hl(z)Qi(z) 

• roots of i*^ 

• roots of H2 

• root of Q l 
O roots of Pq 



-0—9- 



Re z 



Figure 12: The tree obtained in the (2,2) confining branch with a four-valent vertex. It has 
confinement index t = 3. 

• The (2,2) Coulomb Vacua 

The solution of the factorization problem ( |5.4| ) in this branch is parametrized adi 



19 Refer to equation (3.53) in §1. We have set g — > ah, z > z h and A — > A h in that equation. 
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z 2 + (l + a)z + 



(3 + a)(9 + 15a-a 2 + a J 



P G {z)-2A t 



108 

^ 2 (l-a)(3-a) 2 (3 + a) 
108 

2a. 2 M 2a, (3 - a)(9 - 15a - a 2 - a 3 ) 

(z H r + (1 - o)(z H ) + 

; 3 ' V A 3 ; 108 

2a 2 _ (l + a)(3 + a) 2 (3-g) 

1 3 ; 108 



(5.14) 



with a and A satisfying the constraint 



a 5 ( a 2 - 9) 2 = 27 3 A 6 . 



(5.15) 



TVie Rigid Quartic Factorization 



Requiring that the first factor in the either of the two equations in (|5.14j ) has a double 
root leads to the quartic factorization (|5.4j). We get the condition 



a / -25 = 0. 



(5.16) 



For a = 5, the polynomials that solve the equation (|5.7|) are given by 



Q 1 (z)=z-2, H 2 (z) 
280 



2 2 128 

z 2 z H , F 4 

3 27 ' 



z' + 



64 



20 332 

— z H 

3 27 



P«(z) 



8z° + 



9 



4 800 o 560 2 2048 3839488 
z 4 —z 6 H ^ + 



9 



9 



19683 



(5.17) 

Plotting the zeroes of the polynomials lead to the tree in Figure 13. Since the tree is found 
in the Coulomb branch, it has confinement index 1. This can also be checked directly using 
the combinatorial definition: we get t = GCD(3, 4) = 1. For a = —5, we get an equivalent 
tree but reflected about the Re(z) = axis. Note that the solution in (|5.17| ), like the ones 
we have obtained earlier, are polynomials defined over Q. This is why the Galois orbits 
in each case consist of only a single tree. We now discuss a set of trees whose associated 
number field is non-trivial and therefore constitute a larger Galois orbit. 
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P Q {zf - 4A 12 = F A {z) H 2 2 (z) Qi(z) 

• roots of F4 

• roots of H2 

• root of Q x 
O roots of Pq 




Figure 13: The tree obtained in the (2, 2) coulomb branch with a four-valent vertex. 
The Rigid Cubic Factorization 

Requiring that the two factors in the first of the two equations in ( |5.14| ) have a root 
in common leads to the non-trivial condition!^ 



a s - 3a 2 + 3a + 15 = . 



(5.18) 



20 The similar constraint for the second equation does not change the number field and we get 
the same set of trees. 
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Solving for a leads to three solutions 



a 



= (1-2(2)5) 
a(+) = (l + 2*(l-i>/3)) 
k tr(-) = (l + 23(l + i 73)) 



(5.19) 



Substituting these results into the polynomials and plotting their roots lead to the trees 
in Figure 14, 15 and 16 respectively. All of these have confinement index t = 1. 

From the fact that all three trees are obtained from a single polynomial Q5.18 ) irre- 
ducible over Q, it follows that these three trees belong to the same Galois orbit. From 
( |5.19| ), we observe that unlike the earlier solutions which were all defined over Q, the 
number field associated to these trees is non-trivial. 

Pl{z)-^=F 3 (z)Hj{z)Q\{z) 

• roots of F% 

• roots of Hz 
Im z • root of Qi 

O roots of Pq 



-O • 



-0.75 -0.5 -0.25 



0.25 0.5 0.75 \ 1 




Figure 14: Tree for the case cr^ in ( $.1$) . It has t = 1. 



Let us study this example in more detail. The terminology is explained in Appendix 
A. By inspection of ( |5.19|) , we find that and belong to the field K = Q(2^,u) 
where u> = |(1 + (a cube root of unity). It is easy to see that K is also the splitting 
field of the polynomial x 3 — 2, whose associated Galois group Gal(-RT/Q) is the group of 
permutations of three elements, S3, which is non-abelian. 
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P*(z)-4A^ = F 3 (z)H!(z)Ql(z) 

# roots of F3 
roots of H3 

# root of Qi 
O roots of Pq 




Figure 15: Tree for the case ' in (5.19 ). It has t = 1. 



Im z 




P^z)-4A 12 = F 3 (z)H!(z)Q 3 1 (z) 

# roots of F 3 
roots of H3 

# root of Qi 
O roots of P6 



Figure 16: Tree for the case a^ + ' in (5.19). It has t = 1 
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• The (3,2) And (2,3) Vacuum 

If two roots of F^{z) coincide in (|5.1| ) we get the factorization problem 



P£{z) - 4A 12 = F 2 (z)H£(z) = - ^)H£{z) 



(5.20) 



where, in the second equality, we have suitably shifted and scaled z. Setting A to one, 
( |5.20| ) is solved by the polynomials |T! 



P 6 (z)=2T fl (|) and H s {z) = U 5 (~). 

Plotting the zeroes of these polynomials leads to the branchless tree discussed in Section 
2.8. In Figure 17, we show the tree that arises for the particular case of U(6). As explained 
in the earlier more general discussion, this is the singularity at which the M = 1 branches 
meet. This is the only case we consider that has one of the iVj = 0. Therefore this 
can be explained as the intersection of U(l) and U(l) 2 branches. See Appendix C for a 
more complete discussion. From the combinatorial definition of the confinement index, the 
dessin has t = 6, the largest value for the case with N = 6 edges. 

Im z 

o . 3 • roots of F 2 

# roots of Hk 

0.2 [ 

O roots of Pe 

o.i 



-i 



-0.5 



-0.1 
-0.2 
-0.3 
-0.4 



0.5 



Figure 17: The tree obtained at the maximally confining point. It has t = 6 

5.2. Classifying Dessins From Gauge Theory 

So far we have started with the non-rigid factorization problem and tuned the param- 
eters to get isolated singularities where dessins appear. We have seen how the dessins fall 
into different Galois orbits. We now classify them according to the M = 1 branches to 
which they belong using gauge theory order parameters. We summarize all our findings 
from the gauge theory point of view in Figure 18. Given that the confinement index is a 
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Galois invariant, we find each of the three trees with a 4-valent vertex to belong to distinct 
Galois orbits as shown in Figure 18. Similarly, using the confinement index, we find that 
the trivalent trees fall into at least two distinct Galois orbits: the trivalent dessin with 
t = 2 is left invariant under the action of Gal(Q/Q). 




( S+ , S _) = (3,2) 



Figure 18: Summary of the analysis of the U (6) gauge theory and the location of the trees in the 
various branches of the gauge theory moduli space. The dotted lines indicated a coarse-grained 
classification of J\f = 1 branches based on the (s+,s-) values. In this case the branches are 
more hnely distinguished by the confinement index t. Branches meet at the point where another 
monopole becomes massless. 

Note that in the t = 1 and t = 2 branches we have both trivalent and quartic dessins. 
From our discussion about valency lists, it follows that they are in distinct Galois orbits. 
However, in each case, they have the same value of t and of (s+, s_). This is where extra 
gauge theory criteria are needed in order to distinguish these isolated phases. 

Let us concentrate on the pair of dessins with t = 2. A simple way to see that these 
correspond to two different phases is by tracing them back to the problem in U(3) broken 
to U(l) x £7(2) and then "multiplying by 2". The trivalent dessin comes from the unique 
dessin with three edges. Such a special point is where P$(z) — 1 = z 3 (z 3 — 1). This is 
actually known to be a super conformal field theory in the IR (see [|T3| , |2"4]j and references 

43 



therein). On the other hand, the quartic dessin comes from a generic point in the U(3) 
theory and therefore it is a different phase. 

This discussion proves that the two dessins are in different phases. However, we want 
to go further and show that even the theory corresponding to the quartic dessin is a distinct 
phase from its neighbors in the t = 2 branch. Recall that in Section 4.2 we argued that 
there might be extra chiral ring relations that characterize the corresponding special points 
within each branch. Let us see this in detail in this case. 

For the quartic factorization ( |5.4j ) in the branch defined by = s_ = 2, we take 

H 2 (z) = (z-c) 2 , R 2 (z) = z 2 + az + , H 2 (z) = z 2 + 7 z + 5 . (5.21) 



where we have adopted the notations in ( |3.4j ). Using these in the definition of the t r in 
Section 3.2 we find the following chiral ring relation that is satisfied only at the special 
(quartic) point in the Af = 1 branch: 

4 1 4~3 t\ 4+4 1\ £ 3 -6 1 t x t 2 1 3 +4 t\-4t\ A 6 +6 £ h t 2 A 6 -2tg £3 A 6 +t 2 A 12 = . (5.22) 

Note that this relation uses the A 6 term in ( ]3.6|) . In this polynomial, each term has the 
same R-charge, or equivalently the same dimension, and also the same To see this, 
note that A naturally has Q$> = 1. However, in ( p. 22]) we have set the coupling of the 
cubic term in the superpotential to 1. Such a coupling g is dimensionless, has Q$ = —3 
and shows up in (|5.22|) in the combination gA 6 which then has = 3. 

Our new relation ( |5.22| ) is not satisfied at any other point in the t = 2 branch apart 



from the special point under consideration. At any other point in the same branch one 
can show that all £ r 's with r = 0, . . . , 5 are independent. In other words, one can at most 
find a relation similar to ( |5.22|) , which gives t§ in terms of the other six. Assuming that 
our conjectures about the physical order parameters being Galois invariants are correct, 
this concludes our discussion of the trees from the physics point of view as we have, using 
purely gauge theory criteria, managed to classify the dessins into distinct Galois orbits. 

5. 3. U (6) Gauge Theory: A Mathematician's Point Of View 

We now exhibit how a mathematician would tackle the same problem of classifying 
dessins into Galois orbits. We start with a particular valency list and find all solutions 
to the associated polynomial equations using differentiation methods. In the end, one 
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generically finds a polynomial that factors over Q. Each factor corresponds to a different 
Galois orbit. 

For the U(6) gauge theory perturbed by a cubic superpotential, there are three distinct 
valency lists possible for the rigid factorization problem: 

• Consider the branchless tree shown in Figure 17. From the valency list, one gets the 
polynomial equation 

P 2 (z) - 4 = F 2 (z)H 2 5 (z) = (z 2 - A)Hl(z) , (5.23) 



where, in the second equality, we have suitably shifted and scaled z. In Appendix A of [|16 
we have already shown how to obtain the solution to this equation using the differentiating 
trick. The solutions are Chebyshev polynomials. Plotting the roots of the polynomials, we 
get back the tree in Figure 17. 

• The trees shown in Figures 10, 12 and 13 have the same valency list and arise from 
the polynomial equation 

P 6 (z) 2 -4 = F 4 (z) H 2 (z) Qi(z) . (5.24) 

Differentiating ( |5.24[ ) we get 

2P 6 (z)P^z)=H 2 (z) Ql(z)(Fi(z)H 2 (z)Q 1 (z) + 2F 4 (z)H 2 (z)Q 1 (z)+4F 4 (z)H 2 (z)Q , 1 (z)). 

(5.25) 

Since all polynomials involved are monic, it is easy to see that this leads to two equations 



P^(z)=6H 2 (z) Ql(z), 
12P 6 (z) = Fi(z)H 2 (z)Q 1 (z) +2F A {z)H' 2 {z)Q 1 (z) + 4F 4 (z)H 2 (z)Q' 1 (z) . 

After scaling and shifting the z variable, one can write 

H 2 (z)=z 2 -1, Qi(z) = z + qi , 

6 4 

p 6 (z) = z 6 + Pi > F ^ z ) = ^ + E /< z4 ~ l ■ 

i=l i=l 



(5.26) 



(5.27) 



The first equation in ( |5.26| ) leads to linear equations for the pi, which we can easily solve 
to obtain 

p 3 = 2q x {q\ - 3) , p 4 = -9q 2 , p 5 = -6qf . 
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Substituting this in the second of the two equations in ( |5.26|) leads to 

fi = y9i , h = ^(79 ? ? - 25) , f 3 = ^Qi^ql - 55) , 
h = Y^(" 75 " 718 ^i " 35( ?i) > = ^(25 + 276g? + 7q£ 



(5.29) 



such that gi satisfies the equation 

q 1 (q 2 1 -25)(5qj + 3) = 0. (5.30) 



Each inequivalent solution of ( |5.30| ) leads to a dessin. Thus, each dessin has associated to 



it a specific number field |15[]. In ( |5.30| ), there are solutions obtained by an overall sign 
flip: these do not lead to inequivalent trees. 

Plotting the roots of the polynomials for each of the cases qi = {5, 0, *a/§ } respectively 
leads to the trees in Figures 10, 12 and 13. Since each of the values of q± is a Galois orbit 
in itself (up to an overall sign), the three solutions lead to dessins that belong to different 
Galois orbits. Thus, they should have a different set of Galois invariants. That this is so 
can be checked by computing the monodromy groups. We will postpone further analysis 
to the discussion in Section 5.3. 

• The trivalent trees in Figures 11, 14, 15 and 16 all have the same valency list and 
arise from the polynomial equation 

Pi{z) - 4 = F 3 (z) Hi(z) Ql(z) . (5.31) 

Differentiating the equation as before leads to two equations 



12P 6 (z) = Fi(z)H 3 (z)Qi{z) + 2F 3 (z)H' 3 (z)Q 1 (z) + 3F 3 (z)H 3 (z)Q[(z) 
We can choose to parametrize the polynomials as 



(5.32) 



P 6 (z) = z 6 + J2 P* z6 ~ l > H 3 (z) =z 3 + Y,h 
i=l i=i 

3 

Qi(z) = z, F 3 (z)=z 3 + J2h 



(5.33) 



„ z 

i=l 



where we have used the shift symmetry to set the constant coefficient of Q\ to be zero. 
We will not discuss the solution in detail here, as the analysis is similar to the one we did 
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for the quartic factorization. The solution set is parametrized by (hi, hi) that satisfy the 
relation 

24 h\ - 156 h\ h 2 + 450 h\ h\ - 625 h\ = . (5.34) 

We find two branches of solutions : 

a) hi = hi = : This leads to the simple solutions 

P 6 (z) = z 6 + 2z 3 + ^, H 3 (z)=z 3 + 1, f 3 (z)=z 3 + 2. (5.35) 

The tree that corresponds to this solution is shown in the Figure 11. 

b) hi, h 2 7^ : One can use the scaling symmetry to set hi = 1 and there are three 



solutions which are solutions to the cubic equation for h 2 in (|5.34|) . These are given by 



ho 



4°) = -^(-3 - 2 (2)* + 2§ 



4 +) = ^(6 + 2§(l-tV3)-2(2)s(l + iV3)) • ( 5 - 36 ) 
K h ( ~ } = i(6 + 2§(l+zV3) -2(2)i(l-zV3)) 



The three trees associated to h% , h 2 ^ and h 2 + ^ are shown in the Figures 14, 15 and 
16 respectively. Since h 2 , h 2 ^ are solutions to the polynomial equation (|5.34|) which is 
irreducible over Q, the corresponding dessins are part of the same Galois orbit. Moreover, 
the number field is Q(2 1 / 3 ,w) with w 3 = 1 as it should be from our discussion in Section 
5.1. 

For h 2 = h 2 (Figure 14), we present the polynomials that solve ( |5.32|) : 

D/ . 6 6 5 11250 + 7500 2^- 3750 2! 4 1000 + 4500 2s - 3000 2§ 3 

rp. \Z) — Z H Z + Z H Z 

y ' 5 31250 31250 

44 + 30 23 -512s 
H 



31250 

. > , 2 30 + 202^-102§ 2 + 92^-62§ 
H 3 (z) =z 3 + z 2 + — z + 



125 125 

2 2 -15 + 20 2^-10 23 . 8-62^ 



F 3 (z) = z 6 + -z z H — z + 



5 125 125 

(5.37) 

By direct computation, we have therefore classified into Galois orbits the class of trees 
with 6 edges that we considered in Section 5.1. The results of the mathematical analysis of 
the rigid factorizations are summarized in Figure 19. These coincide with the classification 
we obtained from the gauge theory analysis. 
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P£(z)-^ = F 3 (z)Qi(z)Hi(z) 




t = 2 




P b 2 (z)-4A li = F 4 (z)Qt(z)Hi(z) 




t = 3 



V + + 



Figure 19: Summary of the analysis of the U(6) gauge theory and the associated dessins. We 
have also included the results of gauge theory analysis and indicated the confinement index and 
refined valency list of each figure. 



5.4- Using Galois Invariants To Classify Dessins 

In the previous section, we have shown explicitly how dessins are organized into Galois 
orbits. We now attempt to rediscover the classification using the Galois invariants discussed 
in Section 3.1, focusing on the factorizations in ( |5.24| ) and ( |5.31| ). 

• Quartic Case 

From the direct solution of the factorization problem, we see that there are three 
distinct Galois orbits that correspond, respectively, to the three distinct trees in Figures 
10, 12 and 13. These three trees can been partially distinguished by the more refined 
valency list introduced in Section 2.6 for trees. Let us see this in detail. 
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Depending on how one distributes the roots between the two factors (P(z) ± 2), there 
are two distinct possibilities: 

f Hf(z)Qi(z) 

P 6 (z)-2={ (5.38) 
( F 2 (z)Qf(z) 

If we assign negative valences to each of the zeroes of the polynomials appearing on the 
right, then, the first choice singles out the tree in Figure 12 as the only possibility. On the 
other hand, the second possibility is satisfied by the trees in both Figures 10 and 13. We 
assign signs +/— to the vertices in Figure 19 to indicate these two possibilities. 

In order to distinguish the remaining two dessins, we can compute the monodromy 
group of the trees0. Here, we compute the monodromy groups of the trees with the 
refined bi-partite structure, as in Figure 7 of Section 2.6. Taking the difference of the 
two equations in ( |2.14|) we find an auxiliary polynomial equation that leads to a non-clean 
Belyi map, whose pre-images of 1 are the vertices with — valency and whose pre-images of 
are vertices with + valency. The monodromy group for trees is therefore generated by 
cr + /_, which correspond to the permutation of edges around the +/— vertices respectively. 



Figure 


Monodromy g 


;roup 


10 


(S 3 x S3) x 


c 2 


12 


C 2 x S 4 




13 


S 5 




11 


C3XS3 




14,15,16 


S, 





Table 1. Monodromy groups for dessins 
that occur in the pure U (6) gauge theory 
perturbed by a cubic superpotential. 



The groups of all the trees we have encountered in the U(6) example have been 
collected in Table 1. S n is the permutation group of n elements while C n is the cyclic 
group of n elements. The monodromy group turns out to be different for the trees in 
Figures 10 and 13: for Figure 10, we get the monodromy group (S3 x S3) x C2, of order 
72, while for Figure 13, we get the monodromy group £5, of order 120. 

21 All the monodromy groups have been obtained using the GAP software downloaded from 
http: / /www. gap-system.org . 
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• Cubic Case 



The discussion parallels the one for the quartic factorization. The two possibilities of 
distributing the roots 

f Fs(z)QKz) 

P 6 (z) -2=1 _ (5.39) 

(F 2 (z)Hi(z) 

correspond to two distinct refined valency lists that distinguishes the tree in Figure 11 from 
any one of the trees in Figures 14, 15 or 16. In this case, no further invariant is required 
to distinguish them. 

Thus, we find that the classification of dessins into Galois orbits agrees with what 
we obtained in Section 5.2 regarding the classification of isolated phases in gauge theory. 
So far we have only considered dessins that are trees. We now generalize our discussion 
and consider more general dessins. This will highlight some open questions related to the 
phases of gauge theories with flavour. 



6. Gauge Theories With Flavour 

In this section we turn to a discussion of gauge theories with matter. The dessins that 
appear at isolated singularities in the moduli space will no longer be trees. We will mostly 
focus on the curves that were discussed in |TH] , with isolated Argyres-Douglas singularities 



in the moduli space. 

We start with a general discussion of the non-rigid factorization 

P 2 N (z) + B L (z) = Q 2n (z)H 2 N _ n (z) , (6.1) 

where we have exhibited the degrees of the polynomials explicitly. This curve arises from 
a M = 2 U(N) gauge theory with Nf massive flavors broken to M = 1 by a tree level 
superpotential 

WW = Tr W($) + QjrrJ f ($)Q f , (6.2) 
where / and / run over the number of flavors Nf and 

n+l „ l+l _ 

^) = Ef^> m f f (z) = J2™ f f , k z k - 1 - (6-3) 

k=l k=l 
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The Af = 1 vacua, as before, are those for which Q^niz) = W'(z) 2 + f(z) , where f(z) is a 
polynomial such that deg (/) = deg (W'(z))/2 — 1. m(z) is a matrix of polynomials of size 
N f x Nf. 

It turns out that the only information about the superpotential Qym^{<&)Qf which is 
relevant for the curve ( |6.1|) is the polynomial ]25 , 26^ | 



B L (z) = detm(z). (6.4) 

Clearly, plenty of choices of m(z) can lead to the same Bl(z). 

The particular class of dessins we are interested in arise when (z) has only n + 1 
distinct roots. We use our shift and scale symmetry to set Bl(z) to be of the form 

n 

B L (z) = a z mo (z - l) mi Y[(z- Pj ) m > . (6.5) 

i=2 

Two natural ways of obtaining such 5(z)'s are the following: 

• Nf = L = deg Bl(z) and mj(z) a constant diagonal mass matrix with mo masses 
equal to 0, rti\ masses equal to 1, and rrij masses equal to pj. 

• Nf = n + 1 and rrif(z) a diagonal matrix with polynomial entries z m °, (z — l) mi , and 
(z — Pj) mj . 

The former leads to a theory with unbroken M = 2 supersymmetry if there is no 
W(z). Moreover, it has a large flavor symmetry classically. The latter, on the other hand, 
has a very small number of flavors and generically no special flavor symmetry. 

What we now do to obtain the Argyres-Douglas (AD) dessins studied in |T6[ is to 



further tune the masses of the flavors and the parameters of the superpotential to set 

H N - n (z) = Q 2 n(z) R N -3n(z) ■ (6.6) 



This leads to the rigid factorization problem [|T6 



P 2 N (z) + B L (z) = QUz) R 2 N _ 3n (z) . (6.7) 

We will focus on such factorizations in the rest of the section since many explicit solutions 
to this class of factorization problems have already been obtained in [16J. We will exploit 
these solutions to discuss some interesting issues in gauge theory. 
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6.1. U(10) Gauge Theory With Flavour 

Consider the specific case of dessins arising from the factorization problem of the 
second example in Section 2.3: 

P? {z) + az 5 (z - lf(z - tf = Ql(z)Rl(z). (6.8) 

This problem was completely solved in []T6j . There are two inequivalent solutions to (p] 



Borrowing the explicit solutions from []16] one can plot the zeroes of the polynomials as 
before; we have drawn the corresponding dessins in Figure 20 and 21. 

It turns out that the monodromy group distinguishes between the two dessins and 
therefore they belong to different Galois orbits. For the two dessins considered here, the 
relations a\ = 1 and Oq = 1 are satisfied. Let us denote the two monodromy groups 
corresponding to the two dessins by M\ and M% respectively. In order to explicitly write 
down the generators, it is useful to number the half-edges of the dessin, as in the figures 
20 and 21. 

From the definitions, one can check that 

— Mi is generated by 



en = (1, 2)(3, 4) (5, 6)(7, 8) (9, 10)(11, 12)(13, 14)(15, 16)(17, 18)(19, 20) 
(7 = (1, 10, 11)(2, 15, 3)(4, 5)(6, 19, 7)(8, 9)(12, 13)(14, 20, 18)(16, 17) . 

It has order 30720 and has 84 conjugacy classes. 
— M2 is generated by 

cr x = (1, 2) (3, 4)(5, 6) (7, 8)(9, 10)(11, 12)(13, 14)(15, 16)(17, 18)(19, 20) 
(7 = (1, 10, 11)(2, 13, 3)(4, 19, 5)(6, 7) (8, 9)(12, 20, 18)(14, 15)(16, 17) . 



(6.9) 



(6.10) 



It has order 30720 and has 63 conjugacy classes. 

Since Mi ^ M2, the two dessins belong to distinct Galois orbits. We can now ask if 
it is possible, from the gauge theory analysis, to distinguish between them. Note that the 
holomorphic invariants introduced in section 3.2 do not give any information in this case, 
as it is not possible to factorize ( |6.8j ) as P(z) ± f(z), for some polynomial f(z), unlike the 
pure gauge theory case. This is reflected in the fact that from the mathematical point of 
view, it is not possible in general to define a refined valency list when the dessin is not a 
tree. 
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Jf (z) + az 5 (z - l) 5 (z - tf = Ql(z)Rl(z) 



• roots of Q4 
O roots of P 10 

• roots of R 4 
X poles 




Re z 




-1.5 



Figure 20: One of the two dessins arising from the factorization (6.8). Each face is bounded by 
three line segments containing 1 , 2 and 2 edges respectively. The figure to the right is a schematic 
version of the dessin, with the edges numbered to aid the computation of the monodromy group. 

fg,W + - lf(z - tf = Ql(z)Rl(z) 

# roots of Q4 



Im z 
6 



O roots of P 10 
roots of _R 4 
X poles 




Figure 21: The other dessin arising from (6.t). Each face is bounded by 1,1 and 3 edges between 
trivalent vertices. 
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Similarly, although the combinatorial definition of the confinement index still makes 
sense, there is no sense in which confinement is a good order parameter for gauge theories 
with matter. So although t might still be a good Galois invariant, its physical interpretation 



is unclear. Some preliminary analysis of the phases have already been attempted in |27 
We leave a more detailed study of theories with flavor for future work. 



7. Conclusions And Open Questions 

The central theme of this work has been the relation between dessins d'enfants and 
supersymmetric gauge theory. In this section, we discuss some of the results we have 
obtained and list some of the immediate questions and future directions of research that 
have emerged from our analysis. 

We have seen that any clean dessin with iV c edges and Nf + 1 faces can be found 
at an isolated singularity in the moduli space of an M = 2 U(N C ) gauge theory with Nf 
flavours. The particular rigid factorization of the Seiberg-Witten curve that corresponds 
to the isolated singularity is determined by the valency lists of the dessin. Typically there 
are many solutions to this factorization problem. For each such solution, the rigid curve 
determines a rational Belyi map (3(z) whose inverse image of the [0, 1] interval gives a 
dessin D on the sphere. Such an Af = 2 gauge theoretic perspective matches very closely 
the mathematical point of view of obtaining and classifying dessins into Galois orbits by 
direct solution of the polynomial equation. 

From a physics point of view, it is also natural to study what we called the non- 
rigid factorizations (|i~2|). The solutions to these equations are interpreted as the space of 
vacua that preserve Af = 1 supersymmetry when the Af = 2 theory is perturbed by a tree 
level superpotential. One can associate to these vacua disconnected graphs (which are not 
dessins), that come together and join to form a dessin at any of the isolated singularities 
mentioned in the previous paragraph. The difference lies in the fact that these are now 
looked upon as Af = 1 vacua. 

Interestingly, it is this Af = 1 point of view that has a nice counterpart in the more 
refined mathematical approach to the study of dessins, which is to define a complete list of 
Galois invariants that distinguish the dessins that belong to different Galois orbits. Based 
on the examples we have worked out, we have been led to conjecture a relation between the 
mathematical programme of classifying dessins and the physics programme of classifying 
phases of Af = 1 gauge theory. The strongest form of the conjecture states that every 
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Galois invariant is a physical order parameter that distinguishes different phases in the 
gauge theory. 

One example of a Galois invariant that might give new physical information is the 
monodromy group associated to the dessins. In gauge theory, the group usually discussed 
in the context of Seiberg-Witten theory is the S-duality group: for a genus g Seiberg- 
Witten curve, it is an Sp(2g, 7E) group that acts on the Ai and Bi cycles of the Riemann 
surface. These cycles correspond to the edges that go between the filled vertices of the 
dessins (the pre-images of under the Belyi map). However, in general the monodromy 
group of the dessin involves action on half-edges that go between the pre-images of and 
1 and involves the zeroes of P(z), apart from the zeroes of ysw- It is therefore different 
from the S'-duality group, but what exactly the group signifies in the gauge theory is an 
open question. 

In the other direction, one can ask whether the M = 1 gauge theoretic way of finding 
the dessins, first by solving a more general (non-rigid) factorization problem and then 
imposing suitable constraints so that one approaches an isolated singularity in the moduli 
space, is useful from a mathematical point of view. From our proof that the confinement 
index is a Galois invariant, it seems that the answer is yes. More generally we have seen 
that all the order parameters have a simple interpretation in the mathematical literature 
as Galois invariants. We believe that in this direction the correspondence is on much firmer 
ground. 

In Section 6, we discussed dessins that appear in the moduli space of a U(N) gauge 
theory with matter. Much less is known about the possible phases of the corresponding 
Af = 1 vacua. We exhibited two dessins in a Z7 (10) theory with flavour, which, from 
the mathematical point of view belong to distinct Galois orbits as they have different 
monodromy groups. However, from the physics point of view, with the available order 
parameters, it seems that the points where the two dessins appear describe the same phase. 
By this we mean that none of the known order parameters are of any use to distinguish 
between them. If our stronger conjecture that Galois invariants map to order parameters 
is correct, the monodromy group should correspond to an order parameter in physics that 
can distinguish the two special points where the dessins appear. 

Many possible generalizations of our work present themselves. Since the dessins can 
be drawn on any two dimensional topological surface, it should be possible to extend the 
correspondence we have found to dessins drawn on genus g > 1 Riemann surfaces. On the 
gauge theoretic side, it would be very interesting to classify the dessins that appear in the 
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moduli spaces of more general gauge theories, such as SO/Sp gauge groups, quiver gauge 
theories with products of U (N) factors, etc. 

The study of the Af = 1 branches in M uses the Dijkgraaf-Vafa relation between gauge 



theory and matrix models |2q , p9| , |30|| ■ In this relation, two distinct hyperelliptic Riemann 
surfaces emerge |2lJ : the Seiberg- Witten curve and the spectral curve of the matrix model. 
In the examples considered in the text, the spectral curve was equivalent to the reduced 
Seiberg- Witten curve in ( |1.4| ). However this is not true in general as discussed in the 
examples in Appendix C. It is conceivable that one can tune the parameters of the tree 
level superpotential so that the spectral curve develops an isolated singularity, leading to 
a Belyi map. It might be very interesting to study the dessins that arise this way. 

As we have stressed throughout, the main objective of this article was to provide the 
rudiments of a dictionary between the physics of supersymmetric gauge theory and the 
mathematics related to the action of the absolute Galois group on the children's drawings 
of Grothendieck. Much more work needs to be done to fully understand the correspondence 
between these two fascinating fields of study, which we hope will lead to a deeper and fuller 
understanding of the relevant physical and mathematical problems. 
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Appendix A. Field Theory And The Absolute Galois Group 

The main mathematical object that appears in this paper is the absolute Galois group 
Gal(Q/Q). It is of central importance in many areas of mathematics. Since this group is 
not very familiar to physicists, in this appendix we give a short description of the definition 
stated in the text: Gal(Q/Q) is the group of automorphisms of the field Q of algebraic 
numbers that leaves Q fixed. 

Instead of directly studying Gal(Q/Q), which is a group that cannot even be finitely 
generated, we will study the relevant concepts by first reviewing Galois groups of finite 
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order. Just to give an idea of the complexity of Gal(Q/Q) it is nice to mention that math- 
ematicians are considering the possibility that any finite group can arise as a projection of 
Gal(Q/Q); this is the so-called "Inverse Problem of Galois Theory"il. 

Before going into the details about the different elements that enter in the definition 
of Galois groups let us lay down some field theory basis. In this review we will assume 
familiarity with definitions of fields, rings of polynomials and basic group theory (for a 
very basic introduction and more details of the main example in this appendix see |31]] ). 

Let us start by recalling some basic definitions. We say that a field E is an extension 
of a field F, denoted by F < E, if E has a subfield isomorphic to F. Examples of fields and 
extensions are Q < Q(2 1 / 4 , i) < Q. Our first goal is to review the meaning of expressions 
such as Q(2 1 / 4 ,i). 

Given a field F, a natural object to study is the ring of polynomials, F[z], with 
coefficients in F. From now on we will assume that F is a field of characteristic zero or a 
finite field. This is to avoid certain pathologies that can happen otherwise. Of course, our 
final target, which is Q, has characteristic zero. 

An element a G E is called algebraic over F if it is a zero of a polynomial p(z) in F[z]. 
E is called an algebraic extension of F if all its elements are algebraic over F. We will 
only consider algebraic extensions from now on. It turns out that there exists a unique 
monic irreducible (over F) polynomial such that p(a) = (since F is of characteristic zero 
or finite, a can only be a zero of order one. This is an example of one of the possible 
pathologies we have avoided). Since such a p(z) is unique we call it p a (z). The degree of 
p a (z), deg(p a (z)) = n is also called the degree of a in F. 

A special class of algebraic extensions are those with the structure of a vector space 
with basis {1, a, ... , a 11-1 } and coefficients in F. These are called simple extensions and 
are denoted by F(a). In general, if F < E and E is of finite dimension n as a vector space 
over F, we say that E is a finite extension^ of degree \E : F\ = n over F. 

An example is F = Q and a = 2 1 / 4 . Then Q(2 1 / 4 ) is generated by {1, 2 1 / 4 , 2 1 / 2 , 2 3 / 4 } 
since p a (z) = z 4 — 2 has degree n = 4. 

22 The more precise statement is that any simple group S might be the Galois group of a finite 
normal extension of Q and hence there will be a natural restriction map from Gal(Q/Q) onto S. 

23 Quite nicely, in our case (with F = Q), all finite extensions are also simple! This is called 
the primitive element theorem. 
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Note that z 4 -2 is reducible over Q(2 1 / 4 ), in fact, z 4 -2 = (z-2 1 / 4 )(z+2 1 / 4 )(z 2 +2 1 / 2 ). 
The last factor, (z 2 + 2 1 / 2 ), is irreducible of degree 2 in Q(2 1 / 4 ). So if we adjoin a root of 
z 2 + 2 1 / 2 to Q(2 1 / 4 ) then z 4 - 2 splits over this new field. 

The new element we need is 2 1 / 4 z. However, multiplying by 2 -1 / 4 G Q(2 1 / 4 ) we get 
i. Therefore, the new field is (Q(2 1 / 4 ))(i) = Q(2 1 / 4 ,i). The latter notation shows the fact 
that the order in which we adjoint 2 1 / 4 and i to Q is irrelevant. 

We have achieved our first goal: Q(2 1 / 4 , i) is called the splitting field of z 4 — 2. More 
generally, an extension E of Q is a splitting field if there is an irreducible polynomial in 
Q[z] such that E is the smallest field that contains all its roots. 

Now we need to introduce the concept of the algebraic closure of a field F. A field 
K is called algebraically closed if every non-constant polynomial in K[z] has a root in K. 
Such a K is called an algebraic closure of F if K is an algebraic extension of F, and it is 
denoted by F. F is unique up to isomorphisms. 

The next goal is to study automorphisms of fields. Splitting fields are important 
because given any one of them, say E such that F < E < F, any automorphism of F that 
fixes F maps E onto itself and induces an automorphism of E leaving fixed F. Moreover, 
splitting fields are the only ones with this property. The basic automorphisms are quite 
simple: if a and b are roots of the same irreducible polynomial, then the map (f>(a) = b 
with (j)(q) = q if q G F is an automorphism, a and b are called conjugates and (p is a 
conjugation. Automorphisms of E that leave fixed F form a group under composition 
denoted by G(E/F). If E is a splitting field, then G(E/F) is called the Galois group of E 
over F and it is denoted by Gal(E/F). 

Let E be a finite extension of F. The number of isomorphisms of E into F leaving F 
fixed is the index {E : F} of E over F. It turns out that for the fields we consider 

{E : F} = \E : F\ = \G(E/F)\ (A.l) 

where \G(E/F) \ is the order of the group. 

The next step in any algebra book would be to define separable extensions. However, 
over Q, all extensions are separable and we do not need to worry about that. A very 
special role is played by (separable) splitting fields which are then called finite normal 
extensions^. 

24 The parenthesis around "separable" are there to indicate that it can freely be removed. 
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Now we are ready to discuss the Fundamental Theorem of Galois Theory. The theorem 
states that if E is a finite normal extension of F then there is a one to one correspondence 
between intermediate extensions of F and subgroups of Gal(E/F). The correspondence 
is the following: to each extension B of F such that B < E, one associates the largest 
subgroup Gb of Gal(E / F) that leaves B fixed. Moreover, B is a finite normal extension 
of F if and only if Gb is a normal subgroup. In fact, Gal(B/F) is isomorphic to the factor 
(or quotient) group Gal{E / F) /Gb- 

Let us apply this to our example E = Q(2 1 / 4 , i). As discussed above, E is the splitting 
field of z 4 — 2. It has a basis {1, a, a 2 , a 3 , i, ia, ia 2 , ia 3 } where a = 2 1 / 4 . 

Since \E : Q| = 8 we must have |Gal(i?/Q)| = 8. It is a simple exercise to exhibit 
the eight automorphisms of E leaving Q invariant (see section 47.2 of [^T|). Studying the 
composition table one discovers that the group is nonabelian. Moreover, Gal(i?/Q) = D^, 
the dihedral group (the symmetry group of a square). If we denote rotations by kn/2 
(with k = 0,1, 2,3) as p^, mirror images (reflections) as pi and diagonal flips as Si, then 
the identification with automorphisms is collected in the table below (only the action on 
a and i is needed). 





Po 


Pi 


P2 


P3 


Pi 


Si 


P2 


82 


a — > 


a 


i a 


—a 


— i a 


a 


i a 


—a 


—i a 


i — > 


i 


i 


i 


i 


—i 


—i 


—i 


—i 



The lattice of all subgroups of D4 is well known (see section 47.2 of |3~H]). According 
to Galois theory there must be one and only one intermediate extension of Q(2 1 / 4 , i) for 
each subgroup. This gives rise to the lattice of intermediate extensions of Q(2 1 / 4 ,z). Let 
Kh denote the subfield of Q(2 1 / 4 ,i) left fixed by the subgroup H of D4. For example, it 
is easy to check that K^ po P2 t. = Q(v / 2, i)- Note that Q(a/2, i) is also a splitting field and 
hence a finite normal extension. One can easily check that {po, P2} is indeed a normal 
subgroup of -D4! Likewise, consider K{ po fll y = Q(2 1 / 4 ). This is not a splitting field and 
one can check that {po, Pi} is not a normal subgroup of D4. 

Now we can go back to our object of interest: Gal(Q/Q). Here we will follow very 
closely an explanation given in ]2l[] and we will illustrate the main ideas using our example 
of Q(2 1 / 4 ,i). We have already explained the meaning of the algebraic closure of a field. 
Here Q is then the algebraic closure of Q, the field of algebraic numbers. This is clearly 
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a complicated object that can be constructed as the union of all splitting fields over Q, 
which as we know, are finite normal extensions of Q. Let E denote a generic one, then 

Q= |J E (A.2) 
Ee£ 

where S is the set of all such extensions. For each extension E we have the corresponding 
Galois group, Gal(£/Q), of E over Q. 

Consider our favorite example, E = Q(2 1 / 4 , i). Its Galois group over Q is Gal(i?/Q) = 
D 4 . Consider L = Q(y / 2, z). As we saw, L < E. Now, every automorphism of E leaves L 
invariant. This is because L is a splitting field. The Galois group of L over Q is then the 
factor group Gal(L/Q) = Gal(E/Q)/{p ,p 2 }. 

Now there is a natural group epimorphism pe,l '■ Gal(E/Q) — > Gal(L/Q) given by 
the restriction map. That this is an epimorphism, i.e, an onto map, follows from the fact 
that every automorphism of L that fixes Q can be extended to an automorphism of E in 
\E : L\ ways. In our example \E : L\ = 2. Consider the following automorphism of L: 
(V2,i) -> (A -i). This can then be extended to E in two ways a follows: (2 1 / 4 , i) — > 

Consider now Q(z) = K{ po Pl P2 P3 y. This is also a splitting field. We then have the 
following sequence of finite normal extensions Q < Q(i) < Q(\/2 : i) < Q(2 1 / 4 ,z). From 
this sequence we can make an observation that will be very important in the definition of 
Gal(Q/Q): an element (gi,g 2 ) of the cartesian product Gal(Q(z)/Q) x Gal(Q( v / 2, i)/Q) 
can be extended to an element of Gal(Q(2 1 / 4 , i)/Q) if and only if PlMi){92) = 9i- This 
is because if g E Gal(Q(2 1 / 4 , i)/Q) is one of the possible extensions then it has to have a 
consistent action on each of the subfields. 

The absolute Galois group Gal(Q/Q) can now be constructed in a very similar way. It 
is a subgroup of the cartesian product of the Galois groups of all finite normal extensions 
of Q 

Gal(Q/Q) < J] Gal(£/Q) (A.3) 
Eee 

consisting of all elements (gg) G n^ef Gal(i?/Q) (this is an infinite "array" with one entry 
for each E e S) satisfying the constraint that PK 2 ,K 1 (gK 2 ) = 9k x whenever K\ < K 2 - The 
identification of each g e Gal(Q/Q) with the element (gg) implies that gs is the restriction 
of g to E. That this set of restrictions is consistent follows from the condition involving p. 

Finally, the action of g on Q is determined by the action on each finite normal extension 
via gE- This is the basic result we used in section 2.4 where we discussed the action of 
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Gal(Q/Q) on dessins. We said that if g G Gal(Q/Q) then r\ acts on a dessin by acting on 
the coefficients of the Belyi map. In other words, the coefficients of the Belyi map, being 
found as solutions to some set of polynomial equations, belong to a splitting field E and g 
acts via its restriction 

A.l. Glossary Of Terms In The Text 

• An algebraic number is an element a G C that generates a finite extension Q(a) > Q. 

• Q is the field of all algebraic numbers and it is also the algebraic closure of Q. 

• A number field is a finite algebraic extension of Q. 

• A monic polynomial is one whose monomial of highest degree has coefficient 1. 



Appendix B. The Multiplication Map As A Belyi-Extending Map 

It was shown in |Tl],|8| that once a solution to the factorization problem ( J1.2[ ) is known 
for U(N), then it is possible to construct a solution to a similar factorization problem for 
U(tN). Let us first review this construction. 

Consider the factorization problem 

P 2 (z)-4A 2t = F 2 (z)H 2 _ 1 (z). (B.l) 

The solution is given by 



P t (z) = 2AVT t — - , F 2 (z) = z 2 -4 v 2 A\ H^iz) = ^A^U, 



(B.2) 

where rj 2t = 1. Tt(z) and Ut-i(z) are the Chebyshev polynomials of the first and second 
kind respectively, defined by setting z = cos 9 and 

r - I ~. / \ , „. . . 1 dT+ , . sin(t^) ,„ , 

T t (z) = cos(te) U t -!(z) = ~{z) = —V • B - 3 

t az smt) 

This implies that they satisfy the relation 

l-T?(z) = (l-z 2 )Ul 1 (z). (B.4) 
Now suppose we have a solution to the factorization problem 

P 2 N (z) - 4A 2N = F 2n (z)H 2 N _ n (z) . (B.5) 
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Then we can use the solution to ( |B.2| ) to construct a solution to 

P 2 N (z) - 4A 2tN = F 2n (z)H 2 N _ n (z) (B.6) 

as follows: 



( Pn(z) \ 



P tN (z) = 2k™rfT t , F 2n (z) = F 2n (z) 

H tN . n (z) = V ^A N ^H N . n (z) U t _, , A™ = V 2 A 2N . 

This procedure to get exactly t solutions to the U(tN) theory from a given solution of the 



U(N) theory was referred to as the multiplication map (by t) in ||10|| . 

Let us now show that the multiplication by t map can be used to construct new Beyi 
maps from old ones. The simple example of t = 2 has already been discussed in Section 4. 
As discussed in that section, the main point is to use the multiplication map to define the 
new Belyi map as 

AW = l-§# 5 (B.8) 



where P t N(z) is given by the first equation in ( B.7|) and Pn(z) gives rise to a Belyi map 



^) = l-^#- (B.9) 



We would now like to exhibit ( |B.8| ) as the composition of a Belyi-extending map with the 
Belyi map ( p.9| ). For this purpose, let us define the polynomials 

T 2k (z) =M k (z 2 ) and T 2k+1 (z) = z S k {z 2 ) . (B.10) 

This follows from the form of the Chebyshev polynomials ( B.3| ). If we now define the 
Belyi-extending maps 

a 2k (u) = l-Mlil-u) and a 2s+ i(«) = 1 - (1 - u)S 2 k (l - u) , (B.ll) 

one can check that 

t (z) =a t (/3(z)) (B.12) 

is a Belyi map for all integer values of t. As for a 2 , one check that geometrically, the 
multiplication by t replaces a given edge in a dessin by a branchless tree of length t. 
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Appendix C. Complete List Of Trees For £7(6) 

In this appendix we study the problem of realizing all possible connected trees with 
6 edges. Along the way we will mention some interesting points about the structure of 
M = 1 vacua from the matrix model point of view [|22| . 

In section 5, we studied examples a pure U(6) gauge theory deformed by a cubic 
superpotential. However, a cubic superpotential does not allow enough flexibility in the 
non-rigid Seiberg-Witten curve to realize all possible trees with 6 edges. We mentioned in 
the text that if we consider a superpotential W(z) of degree n + 1, then the curve describing 
the Af = 1 vacua of U (N) is given by 

y 2 = P 2 N (z) - AK 2N = (W'(z) 2 + /^( Z ))4_ B ( Z ) (C.l) 

where W'(z) has degree n. 

However, this is true only for vacua with all values of N{ different from zero. In other 
words, if U(N) is classically broken to U(Ni) x . . . x U(N n ). It turns out that if n — s of 
the iVj's are zero then the description of M = 1 vacua is more subtle and it was elucidated 



in [£2] by using matrix model techniques inspired by the Dijkgraaf-Vafa relation |p8|j29| , p0 
The way to treat all cases at once is by introducing another curve, called the matrix model 
curve, y 2 m = W'(z) 2 + f n -i(z). Then, if U(N) is broken classically to J7(iVi) x . . . x U(N a ) 
the factorization problems to be solved are 

y 2 = P 2 (z) - 4A 2N = F 2s (z)H 2 N _ s (z), 
V 2 m = W'(z) 2 + f n _ x {z) = F 2s (z)R 2 n _ s {z). 

This means in particular, that the (3, 2) vacua in Section 5.1 for the U(6) example is 
the intersection of the s = 2 branch with the s = 1 branch, where the low energy gauge 
groups are U(l) 2 and U(l) respectively. This intersection was studied in great generality 
in 
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Note however, that the trivalent and four-valent factorizations in U(6) do not cor- 
respond to points where the s = 2 branch intersects the s = 1 branch. This is because 
both values of Ni are nonzero. More generally, it is clear that isolated singular points of a 
theory with s = p are particular cases of those of a theory with s = n for any p < n, since 
the corresponding branches can intersect. 
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The set of interesting values of n in the case of U(6) is n = 1, . . . , 6.H We now turn 
to the construction of all possible dessins associated to U(6). It turns out that one only 
needs to consider a quartic superpotential, i.e., n = 3. 

Let us prove that n = 3 suffices. Instead of starting with n = 6 and descending all 
the way to n = 3 it is best to use the fact that isolated singular points give rise to all 
possible connected trees with 6 edges. These trees must have seven vertices. All possible 
valency lists are the following (6, 0, 0, 0, 0, 1), (5, 1, 0, 0, 1, 0), (5, 0, 1, 1, 0, 0), (4, 2, 0, 1, 0, 0), 
(4, 1, 2, 0, 0, 0), (3, 3, 1, 0, 0, 0), and (2, 5, 0, 0, 0, 0). Recall that a valence list (m, ...,u 6 ) 
means that there are Uk vertices with valence k. These were obtained by requiring that 
the sum of all it^'s is always seven and that the sum of all it^'s times k is always 12. 

The previous valence lists lead to factorization problems of the Seiberg-Witten curve 
of the form 

F 6 (z)Hf(z), Fs{z)Q\{z)Hl{z), F 5 (z)Ql(z)Hf(z), 
F4{z)Ql{z)H}{z), F A {z)Q\{z)Hl{z), F 3 (z)Q 2 3 (z)H^z) 1 

and F2(z)Ql(z) respectively. 

Now we can prove our claim by simple inspection. All these factorization problems 
are particular points in the space of curves (with n = 3) given by 

y 2 = Pi(z)-4A 12 =F 6 (z)Hl(z). (C.3) 

It is also easy to check that F 6 (z)H$(z), F 5 (z)Ql(z)Hf(z), F 5 (z)Ql(z)H%(z) 
and F A {z)Q\ (z)-£ff (z) cannot possibly be obtained from the cubic superpotential (the 
n = 2 case considered in the text), while all other factorizations - F^^Q^^Hf (z), 
F 3 (z)Q 3 (z)Hi(z) and F2(z)Q\(z) - were considered in Section 5. 

The inherently n = 3 factorizations give rise to only six different trees: 

25 There is a subtlety when n = 6 which is related to the fact that tr $ r in the superpotential 
is not an independent quantity but this will not affect our discussion (see § for more details). 
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X 



F e (z)H°(z) 





F b {z)Q\{z)Hl{z) 



F 5 (z)Ql(z)Ht(z) 



>-< X' 

F 4 {z)Q\{z)Hl{z) 

Figure 22: All possible dessins (and the associated factorizations) that appear in the M = 2 
moduli space of the C7 (6) gauge theory, apart from the ones already discussed in the main body 
of the article. We have omitted the vertices corresponding to the zeroes of P&(z). 
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